NASA Contractor Report 337 6 


TRANSEP: A Program for 
Lift Separated Flow About 



Airfoils 


Leland A. Carlson 


GRANT NSG-1174 
DECEMBER 1980 


FOR EARLY DOMESTIC DISSEMINATION 

Because of its significant early commercial potential, this 
information, which has been developed under a U.S. Gov- 
ernment program, is being disseminated within the United 
States in advance of general publication. This information 
may be duplicated and used by the recipient with the ex- 
press limitation that it not be published. Release of this 
information to other domestic parties by the recipient shall 
be made subject to these limitations. 

Foreign release may be made only with prior NASA ap- 
proval and appropriate export licenses. This legend shall 
be marked on any reproduction of this information in whole 
or in part. 

Review for general release December 31, 1982 





TECH LIBRARY KAFB, NM 



001,5571 


NASA Contractor Report 337 6 


TRANSEP: A Program for High 
Lift Separated Flow About Airfoils 


Leland A. Carlson 
Texas A & At University 
College Station, Texas 


Prepared for 

Langley Research Center 

under Grant NSG-1174 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


1980 



INTRODUCTION 


In the design and analysis of high performance airfoils, aero- 
dynamicists would like not only to compute cruise behavior, but also to 
predict airfoil pressure distributions at high lift, high angle-of-attack 
conditions. Since such situations are characterized by large regions of 
separated flow and are dominated by strong viscous interaction effects, 
inviscid methods are not applicable. Furthermore, subsonic-transonic 
analysis and design computer codes (ref. 1 and 2) typically only include 
the effects of weak viscous interaction and fail to give acceptable 
answers whenever the length of the separated zone exceeds a few percent 
of the airfoil chord. 

However, Barnwell (ref. 3) recently demonstrated that the direct- 
inverse technique could be successfully applied to the low speed high 
lift case. By specifying the separation point, he was able to obtain 
excellent agreement with experimental data by solving the linear equation 
with direct nonlinear boundary conditions (airfoil ordinates specified) 
upstream of separation and inverse boundary conditions (pressure speci- 
fied) downstream of separation. Thus, the question arose — could similar 
results be obtained using the full potential equation coupled with vis- 
cous interaction and letting the separation point and pressure level be 
determined as part of the solution? 

This report will describe a flow model and computer program, called 

TRANSEP, which can be used to determine the flowfield about a low speed 

single element airfoil at high angle-of-attack, high lift conditions with 

massive boundary layer separation. Since the program is a modification 

of the TRANDES code (ref. 2), which includes weak viscous interaction, 

it can also be used for subsonic-transonic airfoil analysis and design. 
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FLOW MODEL AND SOLUTION METHOD 


As indicated previously, the present approach is based upon the 
direct-inverse method as developed in the TRANDES program and its ability 
to use either the displacement surface (airfoil ordinate plus displace- 
ment thickness) or pressure as the airfoil boundary condition. For the 
low speed high angle-of-attack case, the airfoil lower surface only 
experiences weak viscous interaction. The lower surface boundary layer, 
however, frequently has a long laminar run before transitioning to fully 
turbulent flow. Thus, the model needs to include an initial laminar 
boundary layer in its viscous interaction section. On the upper surface 
the boundary layer is also initially laminar, but it quickly becomes 
turbulent in character followed by boundary layer separation and a sepa- 
rated zone which can extend over as much as three-fourths of the airfoil 
surface. Fortunately, however, the separated region at low frees tream 
velocities is characterized by an approximately constant pressure level. 
Consequently, the low speed massive separation problem has been modeled 
as shown on Figure 1 . 

To obtain the inviscid portion of the flowfield, the exact pertur- 
bation potential equation is solved iteratively using a rotated finite 
difference scheme and column relaxation. To include viscous effects, 
the basic approach is to calculate a boundary layer displacement thick- 
ness for weak interaction regions and to use it to correct the location 
of the displacement surface (i.e., airfoil ordinate plus displacement 
thickness, 6*). For the strongly interacting separated zone, the 
pressure is specified and the location of the displacement surface is 
computed by integrating the surface-tangency conditions, with the 



initial conditions specified by the displacement surface slope and 
ordinate at the interface between the two regions. At present, the 
location and slopes of the displacement surface are updated every ten 
relaxation cycles. 

On the lower surface of the airfoil, the flowfield is determined 
using direct boundary conditions (airfoil specified) including the 
effects of weak viscous interaction. On the upper surface, the flow- 
field is also computed directly with viscous interaction up to the 
separation point, which is determined as part of the solution. Down- 
stream of separation, inverse boundary conditions are used, and the 
pressure is assumed to be constant in the separated zone. Studies with 
this model have shown that the separated-zone pressure, which has to 
be computed as part of the solution, must be determined by conditions 
at both the separation point and at the trailing edge and not just on 
conditions in the vicinity of separation. This result is in agreement 
with the conclusion of Gross (ref. 4) that conditions at the downstream 
end of the separation bubble partially determine bubble pressure. 

In the present formulation, the pressure coefficient for the con- 
stant pressure separated zone is computed by the equation 


, = -2( y ITE - ^sep) 

P , Se P X ITE " X sep 


( 1 ) 


as illustrated on Figure 2. Here (j) and <p are the perturbation 

1 i £j S Gp 

potentials at the trailing edge and the turbulent boundary layer sepa- 
ration point, respectively. Note that equation (1) is a small-pertur- 
bation approximation for C ; and while this form probably introduces 

p , sep 

some error into the overall problem, its usage has been found to be 
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simple, accurate, and adequate. 

In principle, the separated region and wake probably should be 
accurately modeled with respect to physical phenomena and details. This 
approach has been taken by other investigators (ref. 4-7) and would typ- 
ically introduce a large number of computational points to model the 
wake region. In the present model, however, the wake region contains 
very few computational points due to the coordinate stretching. Thus, 
it is treated very simply in that it is assumed to be inviscid with a 
constant pressure-trailing edge formed by the upper and lower displace- 
ment surfaces. As will be discussed later, this simplistic approach 
yields results which agree well with experimental data. In fact, numer- 
ical experiments with the present model indicate that the results for 
the pressure distribution and aerodynamic coefficients are primarily 
dependent upon obtaining accurate predictions for the location of the 
separation point and the magnitude of the separated pressure. Appar- 
ently, the details of the wake region are of secondary importance. 
Finally, as shown on Figure 2, the airfoil circulation, T , is modeled 
as the difference in the perturbation potentials at the airfoil trail- 
ing edge, i.e., r = <|> u - <f> L » and not as (4> ITE - <f> L > • 

In the present TRANSEP code the turbulent boundary layer is computed 
using the Nash-Macdonald method (ref. 8) in the same manner as in the 
original TRANDES program (ref. 2). For the laminar portion, the boundary 
layer is computed using a compressible Thwaites method which is a 
slightly modified version of a NASA Langley code originally developed by 
Grumman Aerospace Corporation. These integral methods are efficient and 
reliable and yield excellent predictions for displacement thickness 
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values (ref. 9). Internally, the transition point is determined from 
a Granville type correlation (ref. 10) based upon the difference between 
the local momentum-thickness Reynolds number and the value at the laminar 
instability point combined with the pressure gradient history. Some- 
times, particularly on the upper surface at high angles of attack, 
laminar separation is predicted upstream of the transition point. In 
these cases, the local momentum- thickness Reynolds number is compared 
to an empirical correlation in order to determine if the laminar bubble 
is long or short. If the estimate indicates that the bubble is long, 
the calculation proceeds, but a warning is printed which indicates that 
the results are probably in error. If the bubble is short, its length 
is assumed to be one horizontal Ax grid width (about 3% of chord) and 
the turbulent-flow computation is initiated at the next downstream grid 
point. As will be demonstrated, this very simple model and approach 
usually yield adequate results. 

The calculation procedure used in TRANSEP is the same iterative 
successive-column-relaxation scheme used in the basic TRANCES program 
except that the separation point and separated pressure level are per- 
mitted to vary. A convergence history for a typical case is shown on 
Figure 3. Initially, some oscillation occurs on each grid, but, as can 
be seen, the values quickly converge. Normally four hundred iterative 
cycles are performed on both the medium and fine grids. Also, the 
location of the separation point is held fixed on the fine grid (typ- 
ically 97 x 49 with 130 points on the airfoil) at the x/c location 
computed on the medium grid (49 x 25 with 66 points on the airfoil) . 
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TYPICAL RESULTS 


In the development of the present model and program, the NACA 4412 
and GA(W)-2 airfoils were primarily used as test cases due to the excel- 
lent experimental pressure distributions available for these airfoils 
(ref. 11-12). Some typical results for these cases are shown on Figures 
4 and 5. In both cases, the lower-surface boundary layer remained 
entirely laminar, although results with an all- turbulent lower surface 
boundary layer showed no significant differences. On the upper surface, 
transition with a short separation bubble occurred near the leading 
edge. As can be seen on the figures, comparison with experimental data 
is good with respect to C^, separation point location, separated pressure 
level, and C^. 

Figures 6 and 7 show comparisons with experiment for these air- 
foils of C versus angle of attack at typical Reynolds numbers. Again, 

L 

at least up to C T , the TRANSEP results are sufficiently accurate 
Lmax 

for engineering usage. 

Another aerodynamic coefficient which would be desirable to predict 
theoretically is drag and, in general, its prediction for extensively 
separated flow is difficult. The present investigation has revealed that 
C integration (called CDWAVE in the code) is totally wrong and that the 
Squire-Young approach (ref. 13) evaluated at the trailing edge is con- 
sistently too high. The latter is expected since the boundary layer is 
not accurately computed in the separation zone. On the other hand, the 
present investigation has shown that when the Squire-Young formula is 
evaluated at separation the results are better. 

The best method, however, seems to be empirical. Thus, the present 



TRANSEP code utilizes Squire-Young evaluated at separation but with the 
reference length being the distance to separation instead of the chord. 
Results using this technique are shown on Figure 8. 

Extensive investigations with the present model have shown that 
the theoretical results are primarily dependent upon the accurate 
prediction of the upper surface separation-point location. In the 
Nash-Macdonald integral method, separation is assumed to occur when a 
parameter, SEP, defined as 

SEP = - - (2) 

q ds 

exceeds a specific value, SP. Here 0 is the local momentum thickness, 
and q is the local boundary layer edge velocity. Obviously, agreement 
between the theoretical and experimental values requires an appropriate 
selection of the value for SP in order to obtain the correct separation 
point . 

Theoretically, in the Nash-Macdonald method, separation occurs 
whenever the displacement-momentum thickness ratio, 6 * / 0 , exceeds three. 
At low speeds and low angles of attack, this value usually corresponds 
to an SP of about 0.0055. In ref. 1, however, it was found that at 
transonic speeds, an SP value of 0.0040 was more appropriate. Inter- 
estingly, the present study indicates that accurate results can only be 
obtained at high angles of attack by making SP a function of angle of 
attack, a, Reynolds number, RN, and the low angle of attack lift curve 
slope, C . The present best estimate for this correlation is: 

J-f 
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+ 0.004714337) 


For a £ 15.3° 

-5 

Use: SP = -7.14352 x 10 a + (0.0142857 C 

L 

i 

If: SP > 0.0055, set SP = 0.0055 

For a > 15.3° 

-3 _1 1 

Use: SP = 3.6213784 x 10 + 0.0142875 C T + (-8.4074 x 10 RN 

Li 

. a 

+ 2.1707 x 10 ) (a - 15.3) 0 

where the angle of attack, a, is always in degrees. 


This formula should be viewed strictly as an empirical attempt to force 

the correct separation point behavior into the Nash-Macdonald approach 

as it is formulated in a finite difference grid. Thus, the user should 

anticipate further empirical modification of it if his experience so 

dictates. This alteration should not be viewed as a criticism of the 

Nash-Macdonald method. Nevertheless, it should lead to reasonable 

results for the flow about low-speed airfoils up to C„ , and some- 

Lmax 

times even beyond C T 

Lmax 
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PROGRAM USAGE 


The program is written in FORTRAN IV programming language for use 
on IBM, AMDAHL, CDC, and similar computers. In nonoverlay mode it 
requires less than 256,000 bytes on an IBM machine. Some modifi- 
cation to formats, etc. , may be required to run the program on different 
computer systems. The changes necessary to convert the TRANDES code, 
as given in Appendix C of NASA CR-2821 , to the present TRANSEP version 
are listed in Appendix A. 

INPUT 

The input deck setup is identical to that of TRANDES. However, 
some of the variables have been redefined for the massive separation 
case, and a few new variables have been added to NAMELISTS FINP and 
IINP. The definitions for these new input variables are as follows: 

XIBDLY - The x-location at which upper surface transition is 

assumed to occur. The turbulent boundary layer calcu- 
lation starts at the next grid point. The relationship 
to percent chord is: 

XIBDLY = (% chord - 50.0)/ 100.0 
Used only if ILAM = 0. Ignored if ILAM = 1. 

Default -0.44. 

XLSEP - When IMASS = 0, location at which the trailing edge 
correction procedure begins. Between XLSEP and the 
trailing edge the pressure distribution and the 
displacement surface is modified. Used in this fashion 
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only if ITEUPC and/or ITELWC equal 1. When IMASS = 1, 
set to a number forward of expected separation point. 
Default 0.50. 

XLBDLY - The x-location at which lower surface transition is 

assumed to occur. Same relationship to chord is XIBDLY. 
Used only if ILAM = 0. Ignored if ILAM = 1. 

Default - 0.44. 

RLAX - Relaxation parameter for separated pressure level. 

Sometimes needed to enhance convergence. It is used 
only when ITACT = 1. IMASS = 1. Default 1.0. 

RADUS - Leading edge nose radius, r LE / c * Used only if ILAM =1. 
Default 0.035. 

CLALP - Low angle of attack lift curve slope. Used only if 
IMASS = 1. Default 0.10966. 

ILAM - Boundary layer parameter. If zero, boundary layer is 
all turbulent. If one, boundary layer is laminar- 
turbulent with natural transition. Default 0. 

IPRT1 - Print parameter. If one, perturbation potential values 
printed at the completion of each grid. Default 0. 
IPRT2 - Print parameter. If one, x and y velocities at each 
grid point printed at the completion of each grid. 

Values printed as part of Mach Chart. Default 0. 

IMASS - Massive separation parameter. It should be one for 

massive separation cases and zero for all others. To 


utilize, ITACT must be 1. Default 0. 


OUTPUT 


The program output for each grid in the massive separation mode is 

1 . Heading 

2. Coordinate System. It is printed as I,X(I) followed by J,Y(J). 

3. Mach number, angle of attack and location where direct 
calculation stops 

4. Case number 

5. Case type callouts, i.e., 

• Inviscid analysis case 

• With viscous interaction 

• Massive separation 

6. Input data in namelists FINP and IINP 

7. Airfoil ordinates in direct region 

X - horizontal ordinate, where -0.5 is leading edge and 0.5 
is the trailing edge 
YU - Upper surface ordinate 
YL - Lower surface ordinate 
Upper Slope - Upper surface slope 
Lower Slope - Lower surface slope 

8. Iteration history at ten-cycle intervals 
CIR - Circulation 

DPM _ Maximum <}> correction (absolute value) in the last 

relaxation cycle, with the corresponding (I,J) grid 
location 

NSSP - Number of supersonic points 

DELTAY or DELSTAR - 6* at the trailing edge for massive 



separation, should go to a constant if solution is 


converging. Only changes after first fifty cycles 
on each grid. 

Separation point and separated are printed. 

9. Laminar Boundary Layer Values (every LP cycles) 

X - Horizontal ordinate 

. 11 

CP - Skin Friction Coefficient (the . 1 x 10 value is an 
arbitrary initial value and should be ignored) 

D-STAR - S*/c, non-dimensional boundary layer displacement 


thickness 

D- THETA - 0 /c, non-dimensional momentum thickness 

H - S*/c, shape factor 

RE-THETA - Reynolds number based on 0 

RE-STAR - Reynolds number based on S* 

TM - © 2 Pressure gradient parameter 

10. Laminar Flow Messages 

a. SEPARATION OCCURRED AT X = ; gives x location where 

laminar separation occurred 

b. SHORT BUBBLE FORMED? TRANSITION TO TURBULENT FLOW 

ASSUMED, X = 

c. LONG BUBBLE? LAMINAR STALL MAY OCCUR, X = 

BOUNDARY LAYER CALCULATION WILL BE CONTINUED AS 
TURBULENT BUT ACCURACY OF RESULTS IS ' QUESTIONABLE 

d. BOUNDARY LAYER CONCULATION COMPLETED? NEITHER 
SEPARATION NOR TRANSITION WAS DETECTED 
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11. Turbulent Boundary Layer Values 
X - Horizontal ordinates 

M - Local Mach Number 

DELS - <5*, displacement thickness 

THETA - 0, momentum thickness 

SEP - parameter used to determine if separation occurs 
H - shape factor 

6 *_ d]3 , pressure gradient parameter 

" T w dx 

TAU - t /pu 2 , non-dimensional wall shear force 

w 

CF, D-STAR, D-THETA, H, RE-THETA, RE-STAR, TM (see 9) 

12. Lower Surface Flow Messages 

a. BOUNDARY LAYER CALCULATION COMPLETED? 

b. NEITHER SEPARATION NOR TRANSITION WAS DETECTED 

13. Upper surface flow message gives separation point and C^ 
for separated region 

14. 1st group - x - horizontal ordinates 

YU, YL, location of the surface ordinates in 
the computational plane 
2nd group - x - horizontal ordinates 

YU, YL, slopes of upper and lower surfaces in 
physical plane 

repeat 8) thru 14) until convergence 

15. Final Boundary Layer Results 

YUORIG - Original airfoil upper surface ordinates 
DU - Upper surface displacement thickness 
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SLU - Slope of upper displacement surface 
YLORIG - Original airfoil lower surface ordinates 
DL - Lower surface displacement thickness 
SLL - Slope of lower displacement surface 

16. Pressure distribution on airfoil 

17. Displacement surface ordinates and slopes in physical 
plane 

18. Mach chart, numbers printed are Mach number multiplied 
by 100. Velocities (U,V) may also be printed. 

19. Wave drag coefficient 

20. Plot of Results 

U - Upper surface C 
v P 

L - Lower surface C 

P 

T - Upper displacement surface 

B - Lower displacement surface 

CLCIR - Lift coefficient from circulation 

CL - Lift coefficient from integration of 

C D - CD WAVE + CDF 

CMLE - Moment coefficient about leading edge 
CDF - Drag coefficient by modified Squire-Young 
CMC4 - Moment coefficient about quarter chord 
NOTE: For low speed massive separation, ignore CDWAVE and 

CD. Use CDF only. 

A sample case input and output is provided in Appendix B to aid 
potential users and for program checkout. It should be noted that 
Appendix B only presents results for a medium grid (49 x 25) calculation. 
Fine grid results for the same case are shown on Figure 5. 



CONCLUDING REMARKS 


Some comment should be made concerning the accuracy and appli- 
cability of the present method. Based upon extensive comparisons 
with experimental data, the TRANSEP method should yield accurate 

massive separation aerodynamic coefficients for airfoils at low speeds 

£ 

at Reynolds numbers above 1 x 10 . Such results can be obtained on 
either the medium (default 49 x 25) or fine grids (97 x 49 normally). 

It should also yield reasonable pressure distribution estimates, 
transition points, and separation locations. Due to the finite mesh 
size used in the finite difference formulation, the separation point 
is probably at best only accurate to three percent chord. Consequently, 
at very high lift coefficients the overall program accuracy may be 
somewhat less. Obviously, the current TRANSEP code is primarily 
limited by the ability of its boundary layer techniques to predict 
separation. 

In conclusion, the user should note that TRANSEP retains all the 
features of TRANDES. Thus, it can also be used for the design and 
analysis of high speed and transonic airfoils. 
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APPENDIX A 


Changes to Convert TRANDES Code As 
Defined in NASA CR-2821 to TRANSEP 



°AGE 19, CHANGE LINES 1 THRU 4 TO THE FOLLOWING 
C **************************** T R A N S E P ************************* 
C ********************************************************************** 
c ***** A FORTRAN PROGRAM FOR COMPUTING THE FLOW ABOUT AIRFOILS ******** 
C ***** AT LOW SPEEDS UNDER HIGH LIFT, MASSIVE SEPARATED FLOW ******** 
C ***** CONDITIONS, TRANSEP IS A MODIFIED VERSION OF TRANDES ******** 
C ***** DISCUSSED IN NASA CR-2821. THUS, IT CAN ALSO BE USED ******** 
C ***** FOR TRANSONIC AIRFOIL ANALYSIS AND/OR DESIGN. ******** 

C ***** LELANO A. CARLSON, AEROSPACE ENGINEERING DEPARTMENT ******** 

C ***** TEXAS A&M UNIVERSITY, COLLEGE STATION, TEXAS, 77843 ******** 

C ***** JANUARY 1980 *************************************************** 
C ********************************************************************** 

PAGE 19, INSERT AFTER LINE 5 THE FOLLOWING 
DIMENSION UVEL( 99) • VVEL<99) 

PAGE 19, INSERT AFTER LINE 25 THE FOLLOWING 

2 , LBDLY, XLBOLY 
COMMON/IPT2/I MASS 

C OM MON/ SETCPS/CPSPO , RL AX 
COMMON/ I PT 3/1 SI DE, I LAM ,RADUS 

PAGE 19, INSERT AFTER LINE 27 THE FOLLOWING 

3, XLBOLY ,RLAX,RADUS,CLALP 

PAGE 19, INSERT AFTER LINE 29 THE FOLLOWING 
2, IMASS, ILAM 
2 • IPRTl , IPRT2 

PAGE 19, INSERT AFTER LINE 30 THE FOLLOWING 
CDSPCFO. 0 

PAGE 20, INSERT AFTER LINE 45 THE FOLLOWING 
IPRTl =0 
IPRT2=0 
I L A M=0 
RAOUS=0 .035 

PAGE 20,, INSERT AFTER LINE 54 THE FOLLOWING 
XLBDLY=-,44 

PAGE 20, INSERT AFTER LINE 57 THE FOLLOWING 
RLAX=1.0 
CL ALP— 

PAGE 20, INSERT AFTER LINE 70 THE FOLLOWING 
I MASS=0 

PAGE 21, INSERT AFTER LINE 111 THE FOLLOWING 
IP( IMASS. EQ. I ) PRINT 6005 
5005 FORMAT < 1 H ,3X,*AND MASSIVE SEPARATION*) 

IF( IMASS. EQ.O )GO TO 60 00 

SPP^O.O 

AL°NU=0 .0 

ALPT=ALPDEG 

IF( ALPDEG.LE. 15.3) GO TO 6006 
ALPNU^ALPT- 15 .3 
ALPT=1 5 .3 

6 0 06 SP^-7. 1 4352E-05* ALPT+I 0 . 0 l 42 85 7*CLALP*0 .00 4714337) 
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IF( ALPDEG.LE. 15.31 GO TO 6008 
SPQ=<-8.4074E-1 14RN+2. 1 707E-04 1 * ALPNU 
SP=SP+SPP 
6008 CONTINUE 

IF( SP.GT.0.0055) SP=9. 0055 
IF< SP.LT.O . 004) SP=0.004 

PRINT 6 007 • SP 

6007 FORMAT! 1H • 3X , • SP PARAMETER VALUE IS CALCULATED TO BE •»Ft0.5> 
6000 CONTINUE 

PAGE 22* INSERT AFTER LINE 138 THE FOLLOWING 
IF< IMASS.EQ.l ) GO TO 5001 
GO TO 5002 

5001 IFIMHALF.GT .2 .OR .ITER. GT.50) GO TO 5002 
GO TO 9 

5002 I F ( I l *GT . I TE1 GO TO 9 

PAGE 22. INSERT AFTER LINE 144 THE FOLLOWING 
IF{ IMASS.EQ.O > GO TO 1441 
JB2= J8— 1 

DO 1442 J=J82,JMAX1 

IF( YU( I TS> .GT .EC J) .AND .YU( ITE) .LE.EC J + l > I GO TO 1443 

1442 CONTINUE 

1443 JA=J+1 

IF( JA.LE.JB) JA = JB4- 1 

PHITE = P< I T E * J A— 1 ) + ( P < l TE * J A ) — P I I TE * J A— 1 ) )4<YUC ITE)-E< JA-1 ) )/OE 
1441 CONTINUE 

PAGE 23. INSERT AFTER LINE 192 THE FOLLOWING 
I F ( IMASS.EQ.l ) DELT AY — YU < I TE ) -YUOR I G ( ITE ) 

IF( IMASS.EQ.l ) GO TO 5003 

PAGE 23, REPLACE LINE 194 WITH THE FOLLOWING 

5003 IF( ITER/10410. EQ.ITERI 

PAGE 24, INSERT AFTER LINE 210 THE FOLLOWING 
IF( MHALF.EQ.2 .AND. ITER .LE.50 I GO TO 9006 
IF( ITER/10410 .EQ .ITER. ANO . IMASS . EQ . I ) CALL SHAPE 
IFC IMASS.EQ.O) GO TO 9006 
00 9112 1=11. ITE 

IF(YU( l ) .GE.O.O) GO TO 9112 
YU ( I )=YU( I— 1 ) 

SLU ( I ) =0 .0 
9112 CONTINUE 

PAGE 24, INSERT AFTER LINE 225 THE FOLLOWING 
IF( IPRT1 .EQ .0 ) GO TO 3776 

PAGE 24, REPLACE LINES 226 THRU 232 WITH THE FOLLOWING 
DO 18 JJ=1.JMAX 
J=JMAX+ l-J J 
PRINT 19. J 

19 F0RMATI1H , 'ROW •.IS) 

PRINT 20, (PI I • J) , 1=1 , IMAX) 

20 FORM AT ( l H .10E11.3) 

18 CONTINUE 

PAGE 25, REPLACE LINES 233 ANO 234 WITH THE FOLLOWING 
PRINT 19. J8 

PRINT 20.1 PS I I ) • 1=1 , IMAX) 
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PAGE 25 * INSERT AFTER LINE 234 THE FOLLOWING 
3T76 CONTINUE 

PAGE 25* INSERT AFTER LINE 247 THE FOLLOWING 
IF! IMASS.EQ • 1 > CALL SHAPE 

PAGE 26* INSERT AFTER LINE 276 THE FOLLOWING 
NMACH=2 
JL=2 

JU=JMAX l 
JOUMj= JM AX— 2 

IF{ JDUM .LE.43 ) GO TO 513 

NMACH^l 

JU= JB 

513 CONTINUE 
PRINT 514 

514 FORMAT! //, 38X* * MACH CHART IN COMPUTATIONAL 
1 PLANE-FREE STREAM FROM TOP*./) 

PRINT 5 15*1 MA X 1 

515 FORMAT! 3X* • 1=2. *, 12* • TOP TO BOTTOM*./) 
PRINT 516* JL * JU 

516 FORMAT! 3X. • J= • • I 2* 1 H. . 12. • LEFT TO RIGHT*/) 

PAGE 26. INSERT AFTER LINE 283 THE FOLLOWING 
UVELI J)=U 
V VEL ! J ) —V 

PAGE 26. REPLACE LINE 23S WITH THE FOLLOWING 
PRINT 23.! IONIC! J). J— JL . JU ) 

PAGE 26. INSERT AFTER LINE 285 THE FOLLOWING 
IF! IPRT2 .EQ .0 ) GO TO 500 
PRINT 2878 • !UVEL< J) • J= JL * J U ) 

PRINT 2 878 . ! V VEL ! J ) . J= J L , J U > 

2878 FORMAT! 1H .I5F8.1) i 

PAGE 26. INSERT AFTER LINE 288 THE FOLLOWING 
UVEH J) -0.0 
VVEL! J) =0. 0 

PAGE 26. INSERT AFTER LINE 301 THE FOLLOWING 
UVEL! J ) =U 
VVEL ! J ) —V 

PAGE 27, INSERT AFTER LINE 316 THE FOLLOWING 
UVEL ( J ) — U 
VVEL! J ) =V 

PAGE 2 7 • REPLACE LINE 318 WITH THE FOLLOWING 
PRINT 28.! IONIC! J).J=JL,JU) 

PAGE 27. INSERT AFTER LINE 318 THE FOLLOWING 
IF( IPRT2 *EQ .0 ) GO TO 502 
PRINT 2878. IUVEL! J) . JJ-. JL/ ) 

PRINT 2878 , ! VVEL ( J ) • J— J L ^ J U ) 

PAGE 27, INSERT AFTER LINE 328 THE FOLLOWING 
UVEL! J ) =U 
VVEL! J)=V 

PAGE 27. REPLACE LINE 330 WITH THE FOLLOWING 
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PRINT 28* ( ION ICC J) • J=JL. JU > 

PAGE 27. INSERT AFTER LINE 330 THE FOLLOWING 
IF( IPRT2.EQ.0 ) GO TO 510 
PRINT 2878. <UVEL< J) • JL % JU ) 

PRINT 2 878. < VVEL( J> * OLj JU > 

PAGE 27. INSERT AFTER LINE 332 THE FOLLOWING 
IF< NMACH.EQ.2 > GO TO 519 
NMACH=NMACH+1 
JL = JB 
JU=JMAX 1 
JOUM=JU-JL-42 
IF( J DUM . GT • 0 ) JU=JL+42 
GO TO 513 
519 CONTINUE 

PAGE 29, INSERT AFTER LINE 400 THE FOLLOWING 
I F ( IMASS.EQ.1 • AN0.MHALF.E0.3) MtTER=400 

PAGE 30. INSERT AFTER LINE 443 THE FOLLOWING 
2 .L9DLY, XLBOLY 

PAGE 30. INSERT AFTER LINE 450 THE FOLLOWING 
L3DLY= I LE-1 
218 L9DLY=LBOLYf 1 

IF( X(LBOLY) .LT.XL9DLY) GO TO 218 

PAGE 33, CHANGE LINE 555 TO THE FOLLOWING 
IF( I .LT .LBOLY ) GO TO 37 

PAGE 35. INSERT AFTER LINE 606 THE FOLLOWING 
INTEGER SEPMK 

PAGE 35, INSERT AFTER LINE 620 THE FOLLOWING 
2 • L8DLY . XLBOLY 
C0MM0N/IPT2/I MASS 
COMMON/F IX/MH ALF 
COMMON/SETCPS/CPSPO.RL AX 
COMMON/ IPT3/I SIDE. ILAM .RADUS 

PAGE 35, INSERT AFTER LINE 625 THE FOLLOWING 
IF( MHALF.EQ.2 . AND.ITER ,LE .50 ) LSPOLD-ITE 
I F ( MHALF.GT.2 .AND. I TER .LE.50) L SPOLD=LSPOLD* 2- 1 

PAGE 35, INSERT AFTER LINE 638 THE FOLLOWING 
JB2— JB— 2 

DO 3010 J=JB2,JMAX1 

IF( YU( I TE) .GT .E( J) .AND .YUC ITE) .LE.E( J+l J ) GO TO 3011 

3010 CONTINUE 

3011 JA= J-frT 

I F ( JA.LE.JB) JA=JB+1 

PHITE=P( ITE, JA-l M-(P( I TE. JA )— P< ITE. J A-l) )*(YU( ITE)-EC JA- l »>/0E 

PAGE 36, INSERT AFTER LINE 641 THE FOLLOWING 
5061 CONTINUE 

PAGE 36. INSERT AFTER LINE 650 THE FOLLOWING 
I'M IMASS.EQ. 1 ) LSEP= ITE 

PAGE 36, DELETE LINES 652 THRU 655 
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PAGE 36. REPLACE LINE 663 WITH THE FOLLOWING 
EMI J ) = 1 • E— 0 3 

PAGE 36. INSERT AFTER LINE 673 THE FOLLOWING 
DSS ( ILEI=SLOWERt ILE ) 

I F ( ISIDE.EQ • I >DSS( ILE) = SUPPER ( I LE ) 

PAGE 37, INSERT AFTER LINE 683 THE FOLLOWING 
I F( ISIDE.EQ . 1 ) ISTART=IBDLY 
I F( IS IDE .EQ .2 ) ISTART^LBDLY 
C I LA M— l 

I FC ILAM • FQ # 0 ) GO TO 5060 

CALL TH WAIT (THET A,HH, JE.DSS.OUDS. 1ST ART ) 

IF( IS IDE.EQ.l ) I BDLY=I START 
IF( ISIDE.EQ.2) LBDLY=I START 
IF( ISTART.GE. ITE > THET2=THETA 
IF( ISTART.GE. ITE ) GO TO 202 
5060 CONTINUE 

IF( ITER/LP*LP .EQ .ITER) PRINT 1 , RN 

1 FORMAT! IHO. IOX. • BOUNDARY LAYER ANALYSIS FOR REYNOLDS NO. OF ».E10. 
1 3.//.5X. • X« ,9X. • M« ,8X, • DELS* *4 X , • THETA* . 3X , • SEP • , 

21 OX. # H* ,9X, •PI* . 5X, *TAU* ) 

PAGE 37. INSERT AFTER LINE 691 THE FOLLOWING 
IF< ILAM .EQ. 1 >THET1=THETA 

PAGE 37, INSERT AFTER LINE 704 THE FOLLOWING 
IF! TAU.LT. 1 .E-06 )TAU=1 .E-06 

PAGE 37. INSERT AFTER LINE 707 THE FOLLOWING 
IFISEP.GT. 1 .0 ) SEP^l.O 

PAGE 37, INSERT AFTER LINE 708 THE FOLLOWING 
IF( ISIDE.EQ. 1 .AND.IMASS.EQ.il GO TO 3005 

PAGE 37, INSERT AFTER LINE 709 THE FOLLOWING 
GO TO 3006 

3005 IF < X< J+ l) »LT • XL SEP) SEP=SP 

3006 CONTINUE 

PAGE 38, INSERT AFTER LINE 741 THE FOLLOWING 
IF{ IMASS .EQ .1 ) GO TO 3001 
GO TO 3004 

3001 I F( ISIDE.EQ. 2 )GO TO 3004 
IF( J.NE.LSPOLO) GO TO 3015 
EMSTR=E MILS POLO) 

USTR=UE I LSPOLD I 
M STR^HH 

T STR=TH T/I X (LSPOLD) +0. 50 ) 

3015 CONTINUE 

IF! SEPMK.EQ. 1 ) GO TO 3004 
IF( SEPR.GT.SP) LSEP=J 
IFILSEP .LT. ITE) SEP*K= 1 
3004 CONTINUE 

PAGE 38. INSERT AFTER LINE 755 THE FOLLOWING 
IF( IMASS. EQ.l ) GO TO 3002 

3002 CONTHNUE 

°AGE 38, INSERT AFTER LINE 758 THE FOLLOWING 
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I F( IMASS »EO • 1 > GO TO 440 


PAGE 39. INSERT AFTER LINE 761 THE FOLLOWING 
440 CONTINUE 

PAGE 43. INSERT AFTER LINE 924 THE FOLLOWING 
IF( IMASS.EQ.O) GO TO 3007 
IFC WHALF.EQ *4 ) LSEP=LSPOLO 
IF (LSEP .GT.LSPOLD)LSEP=LSPOLD 
IFILSEP.GE. C ITE-D) GO TO 3008 
I I =LSEP 
111=11-1 
X1=X( I 1 >-0*001 
X 2= 0*5 
J B2 = JB— 2 

00 3012 J= J82 . J M AX 1 

IF( YU(I 1 >.GT.E< J> *ANO. YU< I 1 > .LE.EC J+l>> GO TO 3013 

3012 CONTINUE 

3013 JA=J+1 

IF( JA.LE.JB) JA=JB+1 

PHSEP=P ( II , JA-1 > +<P< l 1 . JA)-PI II . JA-l ))*(YU( I 1 >-E< JA-1 > >/0E 
CPU ( LSE P) = -2.*(PHITE-P HSEP > / 1 XC I TE) - X< LSEP > > 

CPU( LSEP)=CPSPO+RLAX*< CPU (L SEP > -CPSPO > 

CPSPO=CPU< LSEP) 

DO 3009 J=LSEP,ITE 
3009 CPU( J )=CPU< LSEP) 

PRINT 6 1 52 • X f LSEP ) . C PU € LSEP > 

6152 FORMAT ( • • , 2X • • SEPAR AT I ON AT • . F 1 O • 5 * 5X • • SEP AR ATED CP IS •*F10.5) 
LSPOLD-LSEP 
RETURN 

3008 I1=ITE+1 
111=11-1 
X 1=0.5 
X2=l 0000.0 
3007 CONTINUE 

PAGE 43. INSERT AFTER LINE 945 THE FOLLOWING 
2.LBDLY. XLBOLY 


PAGE 52, INSERT AFTER LINE 
IF< D< J) .EQ.O.O) GO TO 

1292 

21 

THE 

FOLLOW! NG 

PAGE 62. INSERT AFTER LINE 
COMMON/ I PT 2/1 MASS 

1647 

THE 

FOLLOWl NG 

PAGE 63, INSERT AFTER LINE 1722 
IF( IMASS.EQ.l ) GO TO 3600 

THE 

FOLLOW! NG 

PAGE 66, INSERT AFTER LINE 
COMMON/IPT2/I MASS 

1805 

THE 

FOLLOWl NG 

PAGE 68, INSERT AFTER LINE 
I TMP1 = I 1 
ITMP1 1=1 1 1 

I F < IMASS.EQ.l >11 =1 TE-M 
IF( IMASS.EQ . 1 >1 1 1 = ITE 

1876 

THE 

FOLLOWl NG 

PAGE 69, INSERT AFTER LINE 

I i=itmpi 

1939 

THE 

FOLLOWING 


I I l = ITMPl 1 
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noon 


PAGE 70* INSERT AFTER LINE 1956 THE FOLLOWING 
COMMON/ IPT2/I MASS 

PAGE 70* INSERT AFTER LINE 1957 THE FOLLOWING 
DIMENSION PTEMP<99> 

I MARK — 0 

IF< IMASS.EQ.O) TO TO 100 
J92=J0- 1 

DO 101 J=JB2*JMAX1 

I F ( YU( I TE) • GT • E ( J ) • AND • YU C l TE ) ■ LE. E < J+l ) ) GO TO 102 

101 CONTINUE 

102 J A— J* 1 

IF( JA.LE. JB ) JA=J8+1 

STE=S4*2 */P I* AT AN< ( 0 .5-X4 )/A2 ) 

QUAN1-STE-S4 

F=PI*A2*0.5*< l .+T AN( PI 24QUAN1 )**2 ) 

F=1 */F 
JAM 1 = J A — 1 
DO 103 J=J B * J AM 1 
I MARK= 1 

C PTEMP( J) =P< ITE. J ) 

P( I TE * J ) “0 • 5/F*C PU ( ID M4*4P< ITE + 1 • J >-P< ITE+2. J > >/< 2.*OS) 
l +{ STE-SI ITE) )*< -2.*P< ITE+l . . J) fP( l TE+2. J )) /lOS*OS) 

P( I TE, J >=P ( ITE* J )/C 1 .5/DS-I STE-SC ITE ) >/( DS*DS > ) 

PI ( J)=P< ITE* J» 

P2 ( J ) =-P ( I TE+ 1 * J )+2 • *° 1 { J ) 

103 CONTINUE 
100 CONTINUE 

PAGE 71* INSERT AFTER LINE 2017 THE FOLLOWING 
IF( IMARK.EQ.O) GO TO 105 
DO 104 J— JB « J AM 1 

104 P( ITE. J >=PTEMPI J) 

105 CONTINUE 

PAGE 72* INSERT AFTER LINE 2039 THE FOLLOWING 
COMMON/ IPT 2/1 MASS 
IF( Il.GE.ITE) RETURN 

PAGE 72. INSERT AFTER LINE 2040 THE FOLLOWING 
!F( IHASS.EQ.l) GO TO 100 


PAGE 72. INSERT AFTER LINE 2041 THE FOLLOWING 
100 CONTINUE 


PAGE 74* INSERT AFTER LINE 2107 THE FOLLOWING 
IF< IMASS.EQ .0 ) GO TO 1000 
RETURN 

1000 CONTINUE 


PAGE 34, INSERT AFTER LINE 2519 THE FOLLOWING 

SUBROUTINE TH WA I T ( THET A • HH . UE .D SS , DUDS. I START) 

C ***** THIS SUBROUTINE IS BASED UPON A NASA LANGLEY ****************** 
C ***** PROGRAM ORIGINALLY DEVELOPED BY THE GRUMMAN ****************** 

C ***** AEROSPACE CORPORATION ****************************************** 
REAL M 

COMMON CPU ( 99 ) ,CPL( 99) , E( 99 ) , DU 1 ( 99 I ,DU2<99) .DL1 ( 99 > .DL2C99) • DC 99 ) 
1 . FF ( 99) .FFP12 (99) ,FFM1 2(99) .FFM1 f 99) *FFM32<99) • 

l PI ( 99 ) ,P2< 99) *PB( 99) ,P( 99*99 ).RS<99 ) ,S< 99) , SUP { 99 ) . SU0C 99 ) .TEMP ( 99 
2 ) . X ( 99) • Y( 99) • YU (99) • YLI99) * SLUC 99) • SLL ( 99 ) . 

3 A 1 , A2 * A I 2. ALP ,C I R.EPS, EPSS , DE. DS* DP • OPM . F • FP 1 2.FM 1 2. FM32 • M * Q I *QI 2. 
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oj iv t- ro >- oj i\j *- 


4W.X1 . X2 . VV J 0 • VV JB 1 * A A J0 1 • A A JB . QQ JB . Q Q JB 1 . UUJ 0 • VV JBP l • QQ JBP1 * A A J8P 1 
5 . Q . QQtUUJBl «P I «P 12* A 22 .All , X4.S4 
COMMON I « ITE* ITE1. ILE. ILEl.Il.il l . ICON* I MAX* I MAXI* I NV . J8 • JA1 • JB 1 • 

1 JMAX. JCON. JMAX1 • NSSP * I W 

COMMON/ NASH /RN. I BOLY . I TACT . YUOR I G< 99 > . YLOR IG < 99 ) • SUPPER! 99 I .SLOWER 
1(99) .DEL <99 ) • DUPOLD < 99 > .DLWOLOC 99) .COF 
COMMON/ DELT AS ITER 

COMMON/ I PTI /XIBDLY.RDEL.RDELFN. RCPB.SP. XSEP. CONV.CPB * XMON . XLSEP. 

I MI TER.LP. I TEUPC . ITELWC.XPC 
2. L0DLY. XLBOLY 
COMMON/ I PT2 /I MASS 
COMMON/FIX/MH ALF 
COMMON/SETCPS/CPSPO.RL AX 
COMMON/ IPT3/I SIDE. ILAM.RADUS 
DIMENSION UEC99) . DSS < 9 9 ) . DUDS < 99 ) 

DIMENSION TWM< 29 ) ,TWL< 29) . TWHC29) .UBAR<500) 

DIMENSION HP( 29) .HPP( 29) . HPPPC 29 ) .PL (29) . PPL < 29 ) . PPPL < 29 ) 

FI ( YY. Y1 .Y2 . Y3.DB)=YY+DB*< Y1+.5*DB*( Y2+DB4Y3/3. ) ) 

F2( A1 . A 2 . DS ) — .5 *DS*( A1 *A21 

DATA TWL/0.500. 0.46 3.0 .40 4.0 .332.0.359. 0.333.0.313. 

0 .291 . 0.268. 0 .2 44.0.220. 0.203.0. 195.0 . 132.0. 168. 

0.153.0. 138.0.1 30.0.1 22 .0.113.0. 1 04 .0 • 095 • 0 • 085. 

0.0 72. 0.056. 0.0 38.0.027.0 .01 5.0. 000/ 

DATA TWM/2.00 • 2 • 07* 2 • 1 8. 2 • 23 • 2 • 28 . 2 • 34 . 2 • 39. 2.44.2.49. 

2.55. 2. 61 . 2.64.2. 67 .2.71 .2.75.2. 31 .2 .87 .2 .90.2. 94 . 2.99. 
3.04.3.09.3. 15.3.22 .3. 30. 3 .39. 3.44.3 .49 .3 .55/ 

DATA TWM/-0. 25. -0.20. -0.14. -0.1 2.-0. 10. -0.080.-0. 064. 

— O . 048. — 0. 032 .—0 . 01 6.0.0.0.008.0.016 .0. 02 4.0.032.0.040. 

0.048.0. 052.0.056 . 0 • 060 • 0 • 064 . 0 • 068 • 0 • 0 72 • 

0.07 6.0.08 0*0 .084 .0 .086 *0.088. 0. 0 90/ 

data mm /o/ 

IBU8=0 

IF(MM.NE.O) GO TO 2 

CALL F I T 2 < 29 .TWM.TWH. HP. HPP.HPPP.3. 3.0 . .0 •) 

CALL FIT2 < 29.TWM.TWL. PL.PPL.PPPL.3.3.0 . .0.) 

MM = l 

RE=RN 
AM ACH=M 
R = 1 .0 
0MEG=1 » 0 
R CV F— 0 • 89 
C RADUS=0.035 

? CONTINUE 
KTRAN= O 
S I NT — O. 

VAL = .2 

UBAR5B= O. 

IF{ ITER/LP4LP .EQ • ITER) PRINT 100 

TRAT s= SQRT< . 5*R A DU S/RE ) * < 1 • + . 2* AMACH**2) **<-.75) 

DEL TA 1 = O. 6 44 74* TRAT 
DELTA2= 0.294 784 TR A T 
H= 2.187 
T M— -0.08695 
RDl = 0. 

RD2= 0. 

CFV=1 .E+10 
M=ILE- 1 
Tl = O. 

R 1 = 1 • 

XSURF=-0.50 

OUOS< ILE)=UE( ILE )/DSS ( ILE ) 
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SLG— DSS < ILE) 

IF( IS IDE.EQ • 1 )GO TO 1001 

00 1002 J= I L E ♦ I T E 

IFC CPU J + l > .LT.CPLC J I > GO TO 1003 

1002 CONTINUE 

1003 IF( J.EQ .ILE)GO TO 1001 
N— J 

XSURF=X <N) 

SLG— DSS C N+ 1 ) 

1001 CONTINUE 

I F ( ITER/LP*LP .EQ. ITER ) PRINT 1 0 1 • XSURF ,CFV * DELTA 1 . DELTA2 *H ,RD2 • 

1 RD1.TM 

5 DELTP1- DELTA1 
D6LTP2- DELTA 2 
U=UE< N+ 1 ) 

UP— DUDS ( N+ 1 ) 

T2 = U**2*( ( l.+.2*AMACH**2*<U**2-l > > >**2.5 

U8AR<N+t)= U 
TMO = TM 

U8 AR5 A— UBAR5B 

QQQ= SQRTC 1 .+ VAL*AMACH**2* ( 1 .-UBARIN+l >**2) > 

FMACH= AMACH*UBAR(N+1 ) /QQQ 
TRAT= 1 .+VAL*FMACH#*2 

UB AR58— UBAR< N+l )**5/TRAT**t 1*5) *R**2 
S I NT— SI NT+F 2 ( UB AR5A * UB AR5B * DSS (N + l ) ) 

TTT= 0. 45*SINT/UBAR ( N+ 1 )**6*TRAT**3/R**2 
TM = -UP4TTT/SQRT (TRAT) 

IF< TM «GT • 0 • 090) GO TO 19 
DO 6 J= 1*28 

IF(TWM( J+l > .GT.TM) GO TO 60 

6 CONTINUE 

60 D3=TM-TWM(J) 

TL = F1 (TWL( J)*PL( J ) . PPL C J ) • PPPL C J > .OBI 

HINC=F1 (TWH( J) • HP( J ) . HPP< J > . HPPP < J ) ,08) 

IF( HINC .LT. 1 .0) HINC- 1.0 
I F ( TL.LT.O. 1 TL= 0. 

H= TPAT* (HI NC + 1 • >-l • 

DELTA2 = SQR T(TTT/RE) *< 1 • + » 2* AM ACH**2 )*♦<— «75l 
DEL T A 1 = H ♦DELTA2 

RD2 = RE*UBAR( N+ II 4DELT A2*QQQ **3 • 

RO 1 — H *RD2 
C F V=2 • * TL/RD2 
DEL ( N + 1 )=0E1_TA1 
N= N+l 
SM=SLG 

SLG=SLG +DSS ( N + 1 ) 

I F ( I TER /LP*LP .EQ. ITER) PRINT lOl.X(N) .CFV, DELTA1 • DELTA2. H. RD2. RD1 • 
1 TM 

14 IF< KTRAN.GT.O ) GO TO 15 

RD2TP= 217.-1 1787 .*TM+ 366762 *TM** 2— 4380 632 #*TM**3 
1 +10453860. *TM**4 

IF( RD2.LT.RD2TR) GO TO 21 
KTR AN= 1 
KT— 1 

1 F ( ITER /LP*LP .EQ.ITER) PRINT 33. X(N) 

33 FORMAT! 3X. • INST ABILITY DETECTED AT X-*.F10.5) 

UNSTM= -TM 
RD2UNS= RD2 
SSHF= 0 • 

SUNST= SM 
GO TO 21 
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15 S SHF— SSHF— F 2(TMG.TM» DS S t N 1 >/CSM-SUNSTJ 

RD2TR= RD2UNS+91 4 • 4-272 50. *SSHF«- 3293 33*SSHF** 2 
I F( RD2. GT.RD2TR 1 GO TO 25 
GO TO 21 
1 9 I BUB— 1 

IF ( ITER/LP*LP .NE • ITER) GO TO 27 
WRITE! 6 • 102 ) X! N) 

% IF! RD2.GT.135.)WRITE!6. 1071 X<N1 

IFIRD2.LT. 1 35. )WRITE (6 • 108IXCN) 

GO TO 27 

21 I F ( N— ITEJ5.28.28 

25 IF! ITER/LP+LP.EQ.ITER) PRI NT 106.XCN) 

27 THET A— OELT A 2 
HH=H 
I START-N 

I F ( I BUB .EQ . 1 ) ISTART=N+1 
RETURN 

23 IF( ITER/LP*LP .EQ. ITER) PRINT 109 
THET A=OELTA 2 
HH = H 

I ST ART = N 
RETURN 

100 FORMAT (12X. 1 X , ,13X. 1 CF • * l 1 X* • D— ST AR • . 9X . * D-THETA 1 »10X."H».12X. 

1 • RE — THET A • ,7X . • RE-STAR* .lOX.'TM'/) 

101 FOPMATI 4X. 8EI 5. 4 ) 

102 FORMAT! / 1 0 X , 1 SEPARATION OCCURRED AT X= *.F10.5> 

C 103 FORMAT ( /10X.4FIN AL LAMINAR BOUNDARY LAYER PROFILES*/! 

C 104 FOPMAT! /10X.4Y* • 1 9X .*VEL. R ATI O* • 1 0 X .* ST RE SS */ ) 

C 105 FORMAT! 3! I OX, FI 0 .5) ) 

106 FORMAT! /10X. • TRANSIT ION OCCURS AT X= '.FIO.S! 

107 FORMAT! 1 0 X * * SHORT BUBBLE FORMED? 1 . 

1 'TRANSITION TO TURBULENT FLOW IS ASSUMED. X= 1 ,F 10 • 5/ ! 

103 FORMAT! 10X,*L0NG BUBBLE? LAMINAR STALL MAY OCCUR. X= 1 ,F10.5. 

1 /I OX. 1 BOUNDARY LAYER CALCULATION WILL BE CONTINUED '. 

2 * A S TURBULENT BUT ACCURACY OF RESULTS IS QU EST l ONABLE 1 / 1 
109 FORMAT! //10X, 'BOUNDARY LAYER CALCULATION COMPLETED? 1 . 

I /10X, 'NEITHER SEPARATION NOR TRANSITION WAS DETECTED 1 /) 

C 110 FORMAT ! //10X, *SP6CIFIED SEPARATION POINT REACHED. X=*.F10.5, 

C 1 / 1 0 X , *NO TURBULENT CALCULATION WILL BE PERFORMED*! 

END 

SUBROUTINE FIT2! N.X.F.FP , FPP ,FPPP ,K 1 .KN.ENOl • ENDN ! 

C IMPLICIT REAL*4 CA-H.O-Z) 

DIMENSION X! I ) .F 1 1 ) , FP! 1 ) .FPP! I ) .FPPP! 1 ) 

NM 1 = N- 1 

IF (N.LT.3) X«N+1) = X!N)4-1.0 

IF (N.LT.3) F(N+1)=F(N) 

DX2-X ( 2 )-X< I ) 

GO TO ! 1,2. 3) ,K1 

1 FP ( 1 ) =0 . 5 

FPPP! 1 )=3. *! ! F( 2 )-F! 1 ) J/DX2-END1 ) /DX2 
GO TO 4 

2 FP!1)=0.0 
FPPP{ 1 ) =END1 
GO TO 4 

3 FP{ I >=-1.0 

FPPP! 1 ) =— OX 2* END 1 

4 DO 5 1=2, NM1 

IP = 1+1 

IM = I - 1 

DX1 = X! I ) — X ! IM ) 

DX2=X ! I P )-X(I) 
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FP ( I )=.5*DX2/(0X l*DX2-.5*OXl*FP< IM ) ) 

5 FPPP ( I)-=s( 6 **(F( IP >-F( I )) /DX2-6**CF( 1>-F| IM 1 )/DXl— DX 1*FPPP( I M )) 
1 ♦FP ( 1 )/OX2 

DX 1 = X ( N )— X ( NMl) 

FPC N J = 0 .0 

GO TO (6*7*8) »KN 

6 FPPP ( M ) =( 6 • *( END N— ( F < M )— F C NMl ) ) /DX1 } /OX I —FPPP ( NW1 ) I /( 2*— FP{ NMl) ) 

GO TO 9 

7 F PP P ( N 1 =ENON 
GO TO 9 

8 FPPP{N) = (ENDN*DXI+FPPP ( NMl ) >/( I • fFP ( NMl > ) 

9 F°P(N)=FPPP(N) 

00 10 II-l.NMl 

1 =N— I I 

IP = 1+1 

DX2=X(IP J-XCIl 

FPP( I >=FPPPt I )-FP( I >*FPPt IP > 

FPPP( I ) = (FPP( IP )-FPP( I))/DX2 
10 FP(I) = (F(IP )-F( 1 ) > /DX2-9X2MFPPI IP ) +2 • #FPP < 1 ) ) /6 • 

FPPP(N) -FPPP C NM 1 ) 

QX1 = X( N)-X< NMl ) 

FP(N)=(F(N)-F(NMl))/OX 1 +DX l ♦ C FPP ( NM l >+2 .♦FPPIN) >/ 6 . 

RETURN 

END 
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APPENDIX B 


Sample Case Input and Output 


SAMPLE CASE INPUT DATA 


NACA 4*12 5 1 ST RUN. MASSIVE SEPERAT I ON 


6FINP M=0 

• 15. ALP* 15 

.3.RN-6 .3+06. XI8DLY= 

-0.47. SP=0 

•00425.XPC 

= 0.4. 


«LAX=0.90, R AOUS— 0» 

0158. CQNV 

= 1 .E-06.XI_SEP=0.0.CLAI_P=0.095. 

&EN0 


SIINP NH ALF — 2. I T ACT 

= 1 • I LAM= 1 . 

I MASS=l . SEND 




O f 

0.0 

0.0 

0.0004 81 

0.008083 

0.001509 

0.010787 

0.002913 

0.013454 

0.00464 1 

0.016057 

0.006628 

0.018567 

0.008806 

0.020960 

0.01991 1 

0.030350 

0.031450 

0 . 037 756 

0.043229 

0. 044071 

0.055172 

0.049642 

0.067237 

0.054647 

0.079399 

0.059195 

0.091640 

0.063353 

0. 1 03947 

0.067173 

0.116310 

0.070688 

0.120721 

0.073928 

0.1411 74 

0.076913 

0.1 53664 

0.079661 

0.166185 

0.082187 

0.170735 

0. 084503 

0.191308 

0.086618 

0.203903 

0.083542 

0.216516 

0.090282 

0.229146 

0.091845 

0.241788 

0.093238 

0.254442 

0.094466 

0.267106 

0.095533 

0.279776 

0.096445 

0.292453 

0.097206 

0.3051 34 

0.097818 

0.31 7816 

0.098287 

0.330504 

0. 090614 

0.3431 89 

0.098803 

0.355074 

0.098858 

0.368556 

0.098780 

0.381235 

0.098572 

0.3939 10 

0.098236 

0.406477 

0.097781 

0.418945 

0.097241 

0.431410 

0 .096625 

0.443872 

0.095932 

0.456331 

0.095166 

0.468785 

0.094327 

0.401235 

0.093417 

0.493680 

0.092439 

0.506121 

0.091392 

0.518556 

0.090278 

0.530986 

0.089100 

0.5434 10 

0.087856 

0.555829 

0.086550 

0.568242 

0.085181 

0.580648 

0.083751 

0.593048 

0.082261 

0.605442 

0.080711 

0 .617829 

0.079102 

0.630209 

0.077435 

0*642581 

0.075711 

0.654947 

0.073930 

0.667305 

0.072092 

0.679656 

0.070199 

0.691999 

0.068250 

0.704334 

0.066246 

0.716661 

0.064188 

0.728900 

0.062075 

0.7412 91 

0.059909 

0.753593 

0.057688 

0.765887 

0.055414 

0.778171 

0.053087 

0.790448 

0.050706 

0.802714 

0.048272 

0.814972 

0.045784 

0.827220 

0.043243 

0.839459 

0.040649 

0.851688 

0.03800 1 

0.863908 

0 • 03S300 

0.87611 7 

0 .032544 

0.888315 

0. 029735 

0.900504 

0.026871 

0.912681 

0.023952 

0.924848 

0.020979 

0.9370 04 

0.017950 

0.949148 

0.014865 

0.961281 

0.011724 

0.973402 

0 • 008526 

0.985511 

0.005271 

0.997607 

0*001958 



0.0 

QA 

1.0 

0.0 

0.0 





0.0 

0.0 

0.000752 

-0.002575 

0.001940 

-0.005006 

0.003538 

-0.007281 

0.005509 

-0.009390 

0.007803 

-0.011325 

0.010359 

-0.013085 

0.013111 

-0.014668 

0.015984 

-0.016079 

0.029504 

-0.020772 

0.042591 

-0.023633 

0.055437 

-0.025555 

0.068118 

-0.026884 

0.080678 

-0.027799 

0.093141 

-0.028408 

0.105525 

-0.028780 

0. 1 1 7843 

-0.028963 

0.1301 05 

-0.028995 

0. 1 42319 

-0.028903 

0.154491 

-0.028707 

0.166626 

-0.028427 

0.1787 30 

-0.028075 

0.190805 

-0.027665 

0 .202857 

-0.027206 

0.21 4887 

-0. 026707 

0.226398 

-0.0261 76 

0.238894 

-0.025620 

0.250876 

-0.025045 

0.262847 

-0.024457 

0.274809 

-0.023859 

0.286763 

-0.023258 

0.298711 

-0.022657 

0.310655 

-0.022060 

0.322596 

-0.021470 

0.334536 

-0.020891 

0.346475 

-0.020325 

0.358416 

-0. 019776 

0.370358 

-0. 019246 

0.382304 

-0.018738 

0.394254 

-0.018254 

0.406312 

-0. 01 7790 

0.41 8469 

-0.017319 

0.430628 

-0.016838 

0.442791 

-0.016349 

0.454958 

-0.015854 

0.4671 29 

-0.015353 

0.479304 

-0.014849 

0.491483 

-0.014343 

0.503668 

-0.013837 

0.515857 

-0.013331 

0.528052 

-0.012028 

0.540253 

-0.012327 

0.552459 

-0.011030 

0.5646 71 

-0.011339 

0.576890 

-0.010853 

0.589115 

-0.010375 

0.601346 

-0.009904 

0.61 3504 

-0.009442 

0.625829 

-0.008989 

0.638081 

-0.008546 

0.65034 1 

-0.008114 

0.662608 

-0. 007692 

0.674882 

-0.007283 

0.687164 

-0.006885 

0.699454 

-0.006499 

0.711751 

-0.006127 

0.724057 

-0.005767 

0.736372 

-0.005421 

0.748694 

-0.005089 

0.761026 

-0.004770 

0.773366 

-0.004465 

0.785715 

-0.004174 

0.798073 

-0.003897 

0.810440 

-0.003635 

0.822816 

-0.003386 

0.835203 

-0.003151 

0.847598 

-0.002930 

0.860004 

-0.002723 

0.872420 

-0.002529 

0 .684846 

-0.002348 

0.897283 

-0.002181 

0.909730 

-0.002026 

0.922188 

-0.001 884 

0.934658 

-0.001753 

0.9471 38 

-0.001634 

0.959631 

-0.001526 

0.972135 

-0.001429 

0 .984651 

-0.001342 

0.997179 

-0.001264 

0.999834 

-0.001249 





0.0 

— 1 • 0 

0.0 

0.0 








N ACA 4412 51 ST RUN • MASSIVE SEPERAT I ON 


X-V GRID SYSTEM 


2 —0.141 OE 01 3 —0.490 OE 00 4 -0.3T06E 00 S -0.2485E 00 6 -0.1247E 00 7 0.0 

3 0.1247E 00 9 0.2435 E 00 10 0.3706E 00 11 0.4900E 00 12 0.I410E 01 

2 -0 .4 26 IE 00 3 -0.1420E 00 4 -0.2303E-07 5 0.1420E 00 6 0.4261E 00 

MACH NO. IS 0.150 ANGLE OF ATTACK IS 15.300 DEGREES 

DIRECT SOLUTION TO 0.50 

CASE NUMBER 100 
INVI5CID ANALYSIS CASE 
WITH VISCOUS INTERACTION 
AND MASSIVE SEPARAT ION 

SP PARAMETER VALUE IS CALCULATED TO BE 0.0049S 
CFINP 

M= 0.14999998 .W = 1.6999998 .XI* 0.50000000 ,X2* 10000.000 .ALP* 0.26703525 . EPS* 0.0 . EPSS 

0.39999998 *X4> 0.48999995 .S4» 2.0000000 .CONV* 0 .99999943E-06 • A 1= 0.24599999 .A2= 0.14999998 .A3* 

699999 . 

RN= 6300000.0 . XI BDLY=-0 . 46999997 .CIR® 0.0 .CDCORR® 0.0 » ROEL* 0.25000000 .RDELFN* 0.12500 

000 • 

$P= 0.4978515 2E-02.XSEP* 0.44000000 ,RCPB= 0.19999999 .CPB- 0.39999998 . XMQN* 0.46999997 .XLSEP* 0.0 

,XPC = 

0.39999998 . XLBDLY=-0 .44 000000 • RLAX* 0.89999998 ,RADUS= 0 . 1 5799999E-01 .CLALP* 0 .94999969E-0 1 

tEND 
Cl INP 

I MAX- 13 • JMAX= 7.IK4SE* 100. INV* 0 .MI TER® 800.NHALF* 2. 1 T ACT* 1 

, ISKP2® 

0.ISKP3* 0.ISKP4® O.ITERP® 0. IREAD* O.LP® 1000.ITEUPC* O.ITE 

LWC= 

0. I MASS® 1*1 LA M= l • 1 PR T 1 = 0.IPRT2* 0 

CENO 

AIRFOIL COORDINATES 

X YU YL UPPER SLOPE LOWER SLOPE 

-0.49000 0.02216 -0.01286 0.98021 -0.64270 

-0.37058 0.07410 -0.02900 0.24890 0.00219 

-0.24852 0.09419 -0.02502 0.09358 0.04866 

-0.12470 0.09869 -0.01903 -0.01674 0.04269 

0.0 0.09192 -0.01399 -0.08421 0.04152 

0.12470 0.07818 -0.00903 -0.13492 0.03668 

0.24852 0.05861 -0.00509 -0.18095 0.02642 

0.37058 0.03380 -0.00256 -0.22597 0.01526 

0.49000 0.00405 -0.00131 -0.27323 0.00631 


1 TERAT I ON 

10 

CIR 

— 

0.23673 

OPM 

a 

0.01788384 

AT 

1 2 

3 

NSSP 

= 

0 

DELTAY 

OR 

DEL STAR 

St 

0.0064 

ITERATION 

20 

CIR 


0 .40439 

DPM 

ss 

0.01465756 

AT 

12 

3 

NSSP 

s 

0 

DELTAY 

OR 

OELSTAR 

a 

0.0064 

ITERATION 

30 

CIR 


0.54029 

DPM 

a 

0.01221043 

AT 

12 

3 

NSSP 

zx 

0 

DELTAY 

OR 

DELSTAR 

2E 

0.0064 

1 TERAT I ON 

40 

CIR 


0 .65274 

DPM 

a 

0.01024419 

AT 

12 

3 

NSSP 

a 

0 

DELTAY 

OR 

OELSTAR 

a 

0.0064 

ITERATION 

50 

CIR 


0 .74703 

DPM 

a 

0.00862175 

AT 

12 

3 

NSSP 

= 

0 

DELTAY 

OR 

OELSTAR 

a 

0.0064 


CP BY CENTRAL DIFFERENCES 
X CPU CPL 


-0.490 

-4.034 

1 .0 05 

-0.371 

-2.849 

0.377 

-0.249 

-1 . 763 

0.330 

-0.125 

-l .232 

0 .327 

0.0 

-0.841 

0.305 

0.125 

-0.555 

0 .291 

0.249 

-0.266 

0.273 

0.371 

0.54 5 

0.244 

0.490 

-0.471 

0.118 



1.410 


0.029 


0.031 


ii n 



OJ 

N3 


-0.49000 

-0.37058 

-0.24352 

-0.12470 

0.0 

0.12470 

0.24352 

0.37058 

0.49000 


0 .0221 6 
0.07410 
0.0941 9 
0.09869 
0.09 192 
0.07816 
0.05861 
0.03380 
0.00405 


YL 

-0.01286 
-0.02900 
-0.02502 
-0.01903 
-0.01399 
-0.00903 
— 0. 00509 
-0.002 56 
-0.00131 


SLU 

0.93021 
0.24890 
0.09358 
-0 .01674 
-0.08421 
-0 . 1 3492 
-0.1 8095 
-0.22597 
-0.27323 


SLL 

-0.64270 
0.00219 
0. 04866 
0.04269 
0. 04152 
0.03668 
0.02642 
0.01526 
0.00631 


MACH CHART IN COMPUT AT IONALTIONAL PLANE-FREE STREAM FROM TOP 

I =>2.12 TOP TO BOTTOM 
J= 2. 6 LEFT TO RIGHT 


15 

15 

15 

15 

16 

7 

2 

0 

28 

25 

12 

1 1 

0 

25 

19 

12 

12 

0 

23 

19 

12 

1 2 

0 

21 

18 

12 

12 

0 

20 

17 

13 

12 

0 

18 

17 

13 

12 

0 

16 

16 

13 

13 

0 

15 

15 

11 

12 

0 

10 

14 

14 

14 

14 

14 

14 


WAVE CO « 0.169714 



CHART 100 


>0.4900 U 

>0.3700 U 

-0.2300 

- 0.1100 

0.0100 

0.1300 

0 • 2500 


TB L 

T B L 

T 8 L 

U T 8 L 

U T B L 

U T B L 

U T B L 


0.3900 


TB L U 


0.4900 


- 4.034 


-3.530 


-3.026 


-2. 522 -2.018 

PRESSURE COEFFICIENT 


CL 


CPSTAR * -29.4074 

1.4907 CO * 0.169714 CMLE = 


- 1.514 


- 1.011 


U 

-0.307 


TB L 
-0.003 


0.501 


1.005 


CLC.IR * 
-0.3039 COF 


l .4041 

0.0 


CMC4 


- 0.0112 


U> 

(jO 


U> 

-O' 


NACA 4412 51 ST RUN. MASSIVE SEPERATION 


X-Y GRID SYSTEM 


2 

-0.3872E 01 

3 

-0.141 OE 

01 

4 

-0. 

a 

-0.3098E 00 

9 

-0.2485E 

00 

10 

-0. 

14 

0.6240E-01 

15 

0. 1247E 

00 

16 

0. 

20 

0.4307E 00 

21 

0. 490 OE 

00 

22 

0. 

2 

— 0.9181E 00 

3 

—0 . 4 26 IE 

00 

4 

-0. 

a 

0.6592E-0I 

9 

0. I 42 OE 

00 

10 

0. 

MACH NO. ISO. 

150 

ANGLE OF , 

ATTACK IS 



DIRECT 

SOLUTION TO 

0 

.50 



CASE NUMBER 100 


1NVISCID ANALYSIS CASE 
WITH VISCOUS INTERACTION 
AND MASSIVE SEPARATION 

SP PARAMETER VALUE IS CALCULATED TO BE 


647 1 £ 

00 

5 

-0.4900E 

00 

6 

-0.4307E 00 

l 86 8E 

00 

1 1 

-0. 1247E 

00 

12 

-0.624OE— 01 

i86ae 

00 

17 

0.2485E 

00 

18 

0.3098E 00 

6471 E 

00 

23 

0. 1410E 

01 

24 

0.3872E 01 

2460E 

00 

5 

-0. 1420E 

00 

6 

-0.6592E-01 

246 OE 

00 

1 1 

0.4261 E 

00 

12 

0.9181E 00 


15.300 DEGREES 


0.00498 


7 -0.3706E 00 
13 0.0 

19 0.3706E 00 

7 -0.9213E-07 


6F INP 

Ms 0.14999998 
0.39999998 
699999 . 

RN= 6300000.0 


.W= 1.6999998 .Xl = 0.50000000 

• X4= 0.48999995 ,S4» 2.0000000 


,X2* 10000.000 .ALPs 0.26703525 

. CQNV* 0 • 9 999994 3 E— 06 » Al« 0.24599999 


• EPS= 0.0 • EPSS 

. A2= 0.14999998 .A3* 


.XI BDLY*-0 .46999997 «CtR= 0.74703497 


,CDCORR= 0.0 


,RDEL= 0. 25000000 * RDELFN“ 0.12500 


000 • 

SP* 0.4978S152E-02,XSEP= 0.44000000 .RCPB= 0.19999999 

,xpc= 

0.39999998 , XL80LY*-0. 44 000000 * RLA X* 0.89999998 

C-END 
61 INP 

IMAX= 25.JMAX* 13.IKASE* 100.INV* 


.CPB= 0.39999998 . XMON= 0.46999997 • XLSEP» 0.0 

• RADUS 3 0 • t 5799999E— 01 *CL ALP* 0 .9499996 9E-0 1 

0 .MI TER* 40O.NHALF* J.I T*Cf* 


•ISKP2= 



0. I SKP3 = 

0 

. ISKP4s 

0. ITERPs 



LWCs 










0. IMASSs 

1 

• I LA M= 

1. IPRTls 




SEND 











AIRFOIL COORDINATES 






X 

YU 

YL 

UPPER SLOPE LOWER SLOPE 




-0.49000 

0.02216 

-0.01 286 

0.98021 

-0.64270 





-0.43067 

0.05546 

-0.02699 

0.38622 

-0.08421 





-0. 37058 

0.07410 

-0. 02 90 0 

0. 24890 

0.00219 





-0.30982 

0.06644 

-0.02769 

0.16175 

0.03592 





-0.24852 

0.09419 

-0.02502 

0 . 09358 

0.04 866 





-0. 1 8677 

0 .098 13 

-0.02193 

0.03544 

0.04963 





-0.12470 

0.09869 

-0.01903 

-0.01674 

0.04269 





-0.06240 

0.09629 

-0.01656 

-0.05562 

0.04025 





0.0 

0.09192 

-0.01 399 

-0.08421 

0.04152 





0.06240 

0.08583 

-0.0114 3 

-0.11 045 

0.04007 





0.12470 

0 . 07818 

-0.00903 

-0.13492 

0.03668 





0. I 8677 

0. 06908 

-0.0069 0 

-0. 1 5826 

0.03197 





0.24852 

0.05861 

-0.00509 

-0.18095 

0.02642 





0.30982 

0 . 04 6 84 

-0.00 36 5 

-0.20333 

0.02069 





0.37058 

0.03380 

-0.00256 

-0.22597 

0.01 526 





0.43067 

0 .01 953 

-0.00179 

-0.24908 

0.01035 





0.49000 

0 .00405 

-0.0013 1 

-0.27323 

0.00631 





ITERATION 

10 CIR = 

0.80321 0PM s 0.00382394 AT 

24 

6 

NSSP 


l TERAT l ON 

20 CIR s 

0.83499 OPM = 0.00253958 AT 

24 

6 

NSSP 

s 

1 TERATION 

30 CIR s 

0.85831 OPM = 0.00207573 AT 

24 

6 

NSSP 

— 

I TERAT I ON 

40 CIR = 

0.87775 DPM s 0.00182289 AT 

23 

6 

NSSP 

— 

ITERATION 

50 CIR s 

0.89464 DPM = 0.00168610 AT 

22 

6 

NSSP 

c 

SEPARATION AT 0 

.37058 

SEPARATED 

CP 1 S 

0.26520 


I TERATION 

60 CIR s 

0.91781 DPM = 0.00202775 AT 

21 

6 

NSSP 


SEPARATION AT Q 

.24852 

SEPARATED 

CP IS 

0 .06477 


ITERATION 

70 CIR s 

0.91161 DPM =.• 0.00149947 AT 

13 

6 

NSSP 


SEPARATION AT 0 

.24852 

SEPARATED 

CP IS 

0.07667 


ITERATION 

80 CIR = 

0.91409 DPM = 0.00119740 AT 

1 2 

6 

NSSP 

s 

SEPARATION AT 0 

. 24652 

SEPARATED 

CP IS 

0 .07555 



0 . I READs O.LPs 

0. IPRT2* 0 


0 DELTAY OR DEL STAR - 
0 DELTAY OR DELSTAR » 
O DELTAY OR DELSTAR * 
0 DELTAY OR DELSTAR = 
0 DELTAY OR DELSTAR * 

0 DELTAY OR DELSTAR = 

0 DELTAY OR DELSTAR = 

0 DELTAY OR OELSTAR = 


1000,1 TEUPCs 


0.0064 

0.0064 

0.0064 

0.0064 

0.0064 

0.0548 

0.1001 

0.1 177 


O.ITE 


U II 



ITERATION 90 CIR 
SEPARATION AT 
ITERATION 100 CIR 
SEPARATION AT 
ITERATION 110 CIR 
SEPARATION AT 
ITERATION 120 CIR 
SEPARATION AT 
ITERATION 130 CIR 
SEPARATION AT 
ITERATION 140 CIR 
SEPARATION AT 
ITERATION 150 CIR 
SEPARATION AT 
ITERATION 160 CIR 
SEPARATION AT 
ITERATION 170 CIR 
SEPARATION AT 
ITERATION 180 CIR 
SEPARATION AT 
ITERATION 190 CIR 
SEPARATION AT 
ITERATION 200 CIR 
SEPARATION AT 
ITERATION 210 CIR 


SEPARATION AT 
ITERATION 220 CIR 
SEPARATION AT 
ITERATION 230 CIR 
SEPARATION AT 
ITERATION 240 CIR 
SEPARATION AT 
ITERATION 250 CIR 
SEPARATION AT 
ITERATION 260 CIR 
SEPARATION AT 
ITERATION 270 CIR 
SEPARATION AT 
ITERATION 280 CIR 
SEPARATION AT 
ITERATION 290 CIR 
SEPARATION AT 
ITERATION 300 CIR 
SEPARATION AT 
ITERATION 310 CIR 
SEPARATION AT 
ITERATION 320 CIR 
SEPARATION AT 
ITERATION 330 CIR 
SEPARATION AT 
ITERATION 340 CIR 
SEPARATION AT 
ITERATION 350 CIR 
SEPARATION AT 
ITERATION 360 CIR 
SEPARATION AT 
ITERATION 370 CIR 
SEPARATION AT 
ITERATION 380 CIR 
SEPARATION AT 
ITERATION 390 CIR 
SEPARATION AT 
ITERATION 400 CIR 
SEPARATION AT 
OJ 
!_n 


■ 0*91772 

0.24852 
s 0.92098 
0.24852 

* 0.922 73 
0.24852 

3 0.92406 

0.24852 
= 0.92549 

0,24852 
= 0.92687 

0.24852 
3 0.92811 

0.24852 
3 0.92923 

0. 24852 
= 0.93024 

0 . 24852 

* 0.93113 
0.24852 

3 0.93194 

0.24852 
3 0.93265 

0.24852 
3 0.93329 


0.24852 
3 0.93386 

0.24852 
s 0 .93438 
0.24852 
= 0.93484 

0.24852 

* 0.93525 
0.24852 

* 0.93561 
0.24852 

= 0.9 3594 

0.248S2 
3 0.93623 

0.24852 
= 0.93649 

0.24B52 
3 0.93672 

0.24852 

* 0.93693 
0.24852 

a 0.93713 
0.24852 
3 0.93729 

0.24652 
3 0.93743 

0.24852 
- 0.93756 

0.24852 

* 0.93768 
0.24852 

= 0.93778 

0.24852 
3 0.93787 

0. 248S2 
= 0.93796 

0.24852 
= 0 .93803 

0.24852 


OPM > 0.00100952 

AT 

13 

6 NSSP 

s 

SEPARATED CP IS 

- 

0.07964 


OPM « 0*00088352 

AT 

13 

6 NSSP 

m 

SEPARATED CP IS 

- 

■ 0 .084 79 


OPM = 0.00078338 

AT 

1 4 

6 NSSP 

s 

SEPARATED CP IS 

- 

0.09480 


DPM = 0.00068170 

AT 

13 

6 NSSP 

3 

SEPARATED CP IS 

— 

0.09909 


DPM 3 0.00059229 

AT 

13 

6 NSSP 

a 

SEPARATED CP IS 

- 

0.10197 


DPM = 0.00051624 

AT 

13 

6 NSSP 

3 

SEPARATED CP I S 

— 

0.10453 


DPM = 0.00045222 

AT 

13 

6 NSSP 

s 

SEPARATED CP IS 

— 

0.10675 


OPM = 0.00039721 

AT 

13 

6 NSSP 

s 

SEPARATED CP IS 

— 

0.10869 


DPM = 0.00034946 

AT 

13 

6 NSSP 

xs 

SEPARATED CP IS 

— 

0.11038 


DPM = 0.00030833 

AT 

13 

6 NSSP 

3 

SEPARATED CP IS 

— 

0.11185 


DPM = 0.00027061 

AT 

13 

6 NSSP 

z 

SEPARATED CP IS 

— 

0.11312 


DPM = 0.0002381 2 

AT 

13 

6 NSSP 

s 

SEPARATED CP IS 

— 

0.11 423 


DPM * 0.00020945 

AT 

14 

6 NSSP 

- 

SEPARATED CP IS 

_ 

0.11519 


OPM = 0.00010471 

AT 

13 

6 NSSP 

X 

SEPARATED CP IS 

— 

0.11 603 


OPM = 0.00016278 

AT 

13 

6 NSSP 

= 

SEPARATED CP IS 

— 

0.11 679 


OPM 3 0 . 000 1 434 I 

AT 

13 

6 NSSP 


SEPARATED CP IS 

— 

0 .1 l 745 


DPM = 0.00012672 

AT 

1 3 

6 NSSP 

= 

SEPARATEO CP IS 

— 

0.11 802 


OPM s 0.00011182 

AT 

14 

6 NSSP 

=5 

SEPARATED CP IS 

— 

0.11853 


OPM = 0.00009835 

AT 

12 

6 NSSP 

X 

SEPARATEO CP I S 

— 

0.11897 


OPM = 0.00008726 

AT 

14 

6 NSSP 

X 

SEPARATED CP IS 

— 

0.11935 


OPM = 0.00007701 

AT 

15 

6 NSSP 

a 

SEPARATED CP IS 

— 

0.11971 


DPM = 0.00006837 

AT 

13 

6 NSSP 

X 

SEPARATED CP I S 

— 

0.12000 


OPM s 0.00006020 

AT 

13 

6 NSSP 

3 

SEPARATED CP IS 

— 

0.12027 


OPM = 0.00005347 

AT 

14 

6 NSSP 

* 

SEPARATED CP IS 

— 

0 .12053 


OPM « 0.00007772 

AT 

24 

5 NSSP 

■ 

SEPARATED CP IS 

— 

0.12075 


OPM = 0.00004196 

AT 

14 

5 NSSP 

* 

SEPARATEO CP IS 

— 

0.12093 


DPM = 0.00003690 

AT 

2 

6 NSSP 

X 

SEPARATED CP IS 

- 

0.12110 


DPM = 0.00003266 

AT 

14 

6 NSSP 

X 

SEPARATED CP IS 

- 

0.12124 


DPM = 0.00002897 

AT 

16 

5 NSSP 

3 

SEPARATED CP IS 

- 

0 .12134 


DPM = 0.00002575 

AT 

16 

6 NSSP 

= 

SEPARATED CP IS 

- 

0.12147 


DPM = 0.00002259 

AT 

13 

6 NSSP 

— 

SEPARATEO CP IS 

- 

0.12158 


DPM = 0.00002056 

AT 

22 

3 NSSP 

= 

SEPARATED CP IS 

- 

0.12168 



0 

DELTAY 

OR 

DELSTAR 

X 

0.1 173 

0 

DEL T AY 

OR 

DELSTAR 

» 

0.1197 

0 

DELTAY 

OR 

OELSTAR 

= 

0.121 a 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1292 

0 

DELTAY 

OR 

DELSTAR 

* 

0.1295 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1293 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1293 

0 

OELTAV 

OR 

DELSTAR 

X 

0.1293 

0 

OELTAY 

OR 

DELSTAR 

- 

0.1294 

0 

DELTAY 

OR 

DELSTAR 

- 

0.1293 

0 

DELTAY 

OR 

DELSTAR 

= 

0 .1293 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1293 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1293 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1293 

0 

DELTAY 

OR 

DELSTAR 

* 

0.1292 

0 

DELTAY 

OR 

DELSTAR 

- 

0.1292 

0 

DELTAY 

OP 

OELSTAR 

X 

0.1292 

0 

DELTAY 

OR 

DELSTAR 

* 

0.1292 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1292 

0 

OELTAY 

OR 

DELSTAR 

E 

0.1 292 

0 

DELTAY 

OR 

DEL STAR 

X 

0.1292 

0 

DELTAY 

OR 

DELSTAR 

X 

0,1292 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1291 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1291 

0 

DELTAY 

OR 

DELSTAR 

- 

0.1291 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1291 

0 

DELTAY 

OR 

DELSTAR 

= 

0.1291 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1291 

0 

DELTAY 

OR 

OELSTAR 

X 

0.1291 

0 

DELTAY 

OR 

DELSTAR 

a 

0.1291 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1291 

0 

DELTAY 

OR 

DELSTAR 

X 

0.1291 



BOUNDARY LAYER ANALYSIS FOR REYNOLOS NUMBER OF 0.630E 07 


X YUORIG DU SLU YLORIG OL SLL 


-0.49000 

0.022 16 

0.0 

0.98021 

-0.0 1286 

0.0 

-0.64270 

-0.43067 

0. 05546 

0.00012 

0.34321 

-0.02699 

0.00005 

-0.08421 

-0.37058 

0.07410 

0.00026 

0.26224 

-0.02900 

0.00012 

0.00219 

-0.30982 

0.08644 

0.00046 

0.16257 

-0.02769 

0.00020 

0.03592 

-0.24852 

0.094 19 

0.00069 

0.09818 

-0.02502 

0.00025 

0.04866 

-0. 1 86 7 7 

0.09813 

0.00096 

0.03992 

-0.02193 

0.00028 

0.04963 

-0. I 2470 

0.09869 

0 .00128 

-0.01 1 15 

-0.01903 

0.00031 

0.04269 

-0.06240 

0. 09629 

0.00168 

-0.04920 

-0.01656 

0.00033 

0.04025 

0.0 

0.09192 

0.00216 

-0.07494 

-0.01 399 

0.00035 

0.04152 

0.06240 

0.08583 

0 .00287 

-0.09601 

-0.01143 

0.00037 

0.04007 

0. 12470 

0. 07818 

0 . 00400 

-0.11273 

-0.00903 

0.00039 

0.03668 

0.18677 

0. 06908 

0 .00S48 

-0.13360 

-0.00690 

0.00040 

0.03197 

0.24852 

0.05861 

0.00709 

-0.10159 

-0.00509 

0.00041 

0.02642 

0.30982 

0.04684 

0.00970 

-0.10519 

-0.00365 

0.00041 

0.02069 

0.37058 

0.03380 

0.01499 

-0.1 0034 

-0.00256 

0.00041 

0.01526 

0.43067 

0.019S3 

0.02307 

-0.07188 

-0.0 0179 

0.00038 

0.010 35 

0.49000 

0.00405 

0.03109 

0. 17174 

-0.00131 

0.00035 

0.00642 


CP BY CENTRAL DIFFERENCES 


X 

CPU 

CP L 



-0.490 

-5.679 

0.738 



-0.431 

-3.487 

0.919 



-0.371 

-2.760 

0.741 



-0.310 

-2.351 

0.660 



-0.249 

-2.000 

0.606 



-0.187 

-1.724 

0.565 



-0.125 

-1.451 

0.525 



-0.062 

-1.188 

0.484 



0.0 

-0.965 

0.449 



0.062 

-0.776 

0.4 19 



0.125 

-0.608 

0.390 



0.187 

-0.401 

0.361 



0.249 

-0.122 

0.329 



0.310 

-0.122 

0.293 



0.371 

-0.122 

0.251 



0.431 

-0.122 

0.201 



0.490 

-0.122 

0.083 



0.647 

0.050 

0.053 



1.410 

0.033 

0.033 



3.872 

0 .022 

0 .022 




YU 

YL 

SLU 

SLL 

-0.49000 

0.02216 

-0.01286 

0.98021 

-0.64270 

-0.43067 

0.05559 

-0.02704 

0.34321 

-0.08421 

-0.37058 

0.07437 

-0.02912 

0.26224 

0.00219 

-0.30982 

0.08690 

-0.02788 

0.16257 

0.03592 

-0.24852 

0.09489 

-0.02526 

0.09818 

0.04866 

-0.1 8677 

0.09909 

-0.02222 

0.03992 

0. 04963 

-0.12470 

0.09997 

-0.01934 

-0.01115 

0. 04269 

-0.06240 

0.09797 

-0.01689 

-0.04920 

0.04025 

0.0 

0.09409 

-0.01434 

-0.07494 

0. 041 52 

0.06240 

0.08871 

-0.011 80 

-0.09601 

0.04007 

0.1 2470 

0.08220 

-0. 00942 

-0.11 273 

0. 03668 

0.1 8677 

0.0746 1 

-0.00730 

-0.1 3360 

0.03197 

0.24852 

0.06720 

—0 .00550 

-0.10159 

0.02642 

0.30982 

0.06039 

-0.004 06 

-0.10519 

0.02069 

0.37058 

0.05417 

-0.00296 

-0.10034 

0.01526 

0.43067 

0.04900 

-0.00217 

-0.07188 

0.01035 

0.49000 

0.05221 

-0.001 66 

0.17174 

0.00642 



I =2 .24 TOP TO BOTTOM 


MACH CHART IN COMPUTAT IQNALTIONAL 


PLANE-FREE STREAM FROM TOP 


J= 2.12 LEFT TO RIGHT 


15 

15 

15 

15 

15 

IS 

IS 

15 

15 

15 

IS 

1A 

14 

14 

14 

15 

15 

15 

15 

IS 

16 

16 

12 

12 

13 

13 

14 

17 

19 

19 

18 

18 

17 

1 1 

9 

7 

5 

4 

0 

33 

28 

26 

22 

19 

12 

1 1 

9 

B 

S 

0 

30 

26 

23 

20 

17 

12 

1 1 

9 

a 

7 

0 

0 

26 

23 

20 

17 

12 

11 

10 

9 

8 

0 

0 

25 

23 

20 

1 7 

12 

11 

10 

9 

9 

0 

0 

24 

22 

20 

1 7 

12 

1 1 

10 

to 

9 

0 

0 

24 

22 

20 

17 

12 

11 

10 

10 

10 

0 

0 

22 

21 

19 

1 7 

12 

It 

11 

10 

10 

0 

0 

21 

20 

19 

17 

12 

11 

11 

11 

11 

0 

0 

20 

20 

16 

17 

13 

12 

11 

11 

1 1 

0 

0 

19 

19 

1 8 

17 

13 

12 

11 

1 1 

1 1 

0 

0 

18 

i a 

18 

17 

13 

12 

12 

12 

12 

0 

0 

17 

1 7 

17 

17 

13 

12 

12 

12 

12 

0 

0 

17 

17 

17 

16 

13 

12 

12 

12 

12 

0 

16 

16 

16 

16 

16 

13 

12 

12 

13 

13 

0 

15 

15 

16 

16 

16 

13 

13 

13 

13 

13 

0 

14 

15 

15 

16 

16 

12 

12 

12 

13 

13 

0 

13 

14 

16 

16 

17 

13 

13 

13 

14 

14 

14 

14 

14 

IS 

15 

16 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

15 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 







WAVE 

CO 

a 

0< 


207096 




CHART 100 


-0.4000 

U 




T8 

L 

-0.4300 


U 



TB 

L 

-0.3700 


U 



T B 

L 

-0.2900 



U 


T B 

L 

-0.2300 




U 

T B 

L 

-0.1700 




U 

T B 

L 

-0.1100 




U 

T B 

L 

-0 .0500 




U 

T B 

L 

0.0100 




U 

T B 

L 

0.0700 




U 

T B 

L 

0.1300 





U T B 

L 

0.1900 





U TB 

L 

0.2500 





UTB 

L 

0.3100 





UTB 

L 

0.3900 





U B L 


0.4500 





U B L 


0.4900 

-5.579 -4.929 

-4.280 -3.630 -2.980 

PRESSURE COEFFICIENT 

-2.330 

-1.680 -1.030 

U B L 

-0.381 0. 

269 0.919 


CL = 1 

CPSTAR = -29.4074 

.7893 CD » 0.221758 CMLE » 

CLCIR 

-0.5022 

- 1.8761 

CDF * 0.014662 CMC4 

« -0.0949 



NACA 4412 51 ST RUN* MASSIVE S6PERAT I ON 


X-Y GRID SYSTEM 


2 

-0.8027E 

01 

3 

-0.3872E 

01 

4 

-0.2276E 

01 

5 

-0.1410E 

01 

6 

-0.91 15E 

00 

7 

—0 .64 71 E 

00 

a 

-0.531 7E 

00 

9 

-0.4900E 

00 

1 0 

— 0 • 460 4 E 

00 

1 1 

-0.4307E 

00 

12 

— 0.4007E 

00 

13 

-0.3706E 

00 

14 

-0.3403E 

00 

15 

-0.3098E 

00 

16 

-0.2792E 

00 

17 

-0.24SSE 

00 

18 

-0.2177E 

00 

19 

-0.186SE 

00 

20 

-0.1558E 

00 

21 

-0. 1247E 

00 

22 

-0.9357E- 

-01 

23 

-0.6240E 

-01 

24 

-0.3121E-01 

25 

0.0 


26 

0.312 1E- 

01 

27 

0.624 OE- 

01 

28 

0. 9357E- 

01 

29 

0. 1247E 

00 

30 

0.1558E 

00 

31 

0.1868E 

00 

32 

0.2177E 

00 

33 

0.2485E 

00 

34 

0. 2792E 

00 

35 

0. 3098E 

00 

36 

0.3403E 

00 

37 

0.3706E 

00 

38 

0.4007E 

00 

39 

0.4307E 

00 

40 

0. 4604E 

00 

41 

0 .4900E 

00 

42 

0.5317E 

00 

43 

0.6471E 

00 

44 

0.91 1 5E 

00 

45 

0.141 OE 

01 

46 

0. 2276E 

01 

47 

0.3872E 

01 

48 

0.8027E 

01 




2 

-0.1 869E 

01 

3 

-0.918IE 

00 

4 

-0.5939E 

00 

5 

-0.4261E 

00 

6 

-0.3206E 

00 

7 

-0 .2460E 

00 

a 

-0.1888E 

00 

9 

—0. 1 42 OE 

00 

10 

-0. 1019E 

00 

1 1 

-0.6592E-01 

12 

-0.3239E- 

-01 

13 

—0*921 3E—07 

14 

0.3239E- 

01 

IS 

0.6592E- 

01 

16 

0. 1019E 

00 

17 

0. 1420E 

00 

18 

0. 188BE 

00 

19 

0.2460E 

00 

20 

0.3206E 

00 

21 

0.4 26 IE 

00 

22 

0.5939E 

00 

23 

0 * 91 81 E 

00 

24 

0. I869E 

01 





MACH NO. IS 0.150 ANGLE OF ATTACK IS 15.300 DEGREES 


DIRECT SOLUTION TO 0.2S 

CASE NUMBER 100 
INVISCID ANALYSIS CASE 
M1TH VISCOUS INTERACTION 
AND MASSIVE SEPARATION 

SP PARAMETER VALUE IS CALCULATED TO BE 0.00498 
tF INP 

Ms 0.14999998 ,W = 1.6999998 .XI* 0.24751765 .X2* 0.50000000 .ALP* 0.26703525 .EPS* 0.0 .EPSS 

0.39999998 ,X4* 0.48999995 «S4* 2.0000000 .CONV* 0 .99999943E— 06 • A 1* 0.24S99999 .A2* 0.14999 998 ,A3= 

699999 • 


RN= 6300000.0 . 

XIBDLYs-0 • 

46999997 

•C1R* 0.93803447 

•CDCORR* 0.0 

• RDEL* 0.25000000 • RDBLFN*' 0*12500 

000 . 







SP= 0. 49785I52E-02. 

XSEP- 0.44000000 . 

RCPBs 0*19999999 

• CP6* 0.39999998 

« XMON= 0.46999997 *XLSEP* 0.0 

.XPC = 







0.39999998 , XLBDLYs-0 . 44 000000 .RLAX= 0.69999998 

•RADUSa 0.15799999E- 

01.CLALP* 0.94999969E-01 

TEND 







61 INP 







1 M AX= 

49 • JMAX* 

2S.IKASE* 

lOO.INVs 

0. MITER* 

400. NHALF* 2,1 TACT* 1 

. 1 SKP2— 








0.1 SKP3= 

0 

,ISKP4= 

0. ITERP= 

0. IREADs 

O.LPs 1000.ITEUPC* 0.1TE 

LWC = 








0*1 MASS* 

1 

* I LA M* 

1 . 1PRT 1 = 

0. IPRT2* 

0 

SEND 









AIRFOIL COORDINATES 




X 

YU 

VL 

UPPER SLOPE LOWER SLOPE 



-0.49000 

0.02216 

-0.01286 

0.98021 

-0.64270 



-0.46044 

0.03887 

-0.02307 

0.66171 

-0.1 9498 



-0.43067 

0 .05359 

-0.02704 

0.3A321 

-0.08421 



-0.40072 

0.06498 

-0.02808 

0.30272 

-0.04101 



-0.37058 

0.07437 

-0.0291 2 

0.26224 

0,00219 



-0.34028 

0.08064 

-0.02850 

0.21240 

0.01905 



-0.30982 

0.08690 

-0.02788 

0. 16257 

0.03592 



-0.27923 

0.09089 

-0.02657 

0.13038 

0.04229 



-0.24852 

0.09489 

-0.02526 

0.09818 

0.04866 



-0*21 769 

0.09699 

-0.02374 

0.06905 

0.04914 



-0.1 9677 

0 .09909 

-0.02222 

0.03992 

0.04963 



-0.1 5577 

0.09953 

-0.02078 

0.01439 

0.04616 



-0. 1 2470 

0.09997 

-0.01 934 

-0.01115 

0.04269 



-0.09357 

0.09897 

-0.01812 

-0.0301 7 

0.04147 



-0.06240 

0.09797 

-0.01689 

-0.04920 

0,04025 



-0.03121 

0.09603 

-0.01562 

-0.06207 

0.04069 



0.0 

0.09409 

-0.01 434 

-0.07494 

0.04 152 



0.03121 

0.09140 

-0.01 30 7 

-0.08548 

0.04079 



0.06240 

0.08871 

-0.01 130 

-0.09601 

0.04007 



0.09357 

0.08546 

-0,0106 1 

-0.10437 

0.03837 



w 0.12470 

0.08220 

-0.00942 

-0.1 1273 

0,03668 




U M 


0.15577 

0.07841 

-0.00 83 6 

-0.12317 

0.03433 





0.18677 

0.0746! 

-0.00730 

-0.13360 

0.03197 





0.21769 

0.07090 

-0.00640 

-0.1 1759 

0.02920 





0.24852 

0.06720 

-0.00550 

-0.10159 

0.02642 





0.27923 

0. 06380 

-0.00478 

-0.10339 

0.02356 





0.30982 

0.06039 

-0.00406 

— 0 • 1 051 9 

0.02069 





0.34028 

0.05728 

-0.00351 

-0.10276 

0.01798 





0.3705B 

0.05417 

-0.00296 

-0.10034 

0.01526 





0.40072 

0.051 68 

-0.00257 

-0.08611 

0.01280 





0.43067 

0 .04900 

-0.002! 7 

-0.07188 

0.0 1 035 





0.46044 

0.05061 

-0.00192 

0.04993 

0.00830 





0.49000 

0.05221 

-0.00166 

0.17174 

0.00642 





ITERATION 

10 CIR a 

0.93756 

OPM a 0.00099355 AT 

38 

2 

NSSP 

X 

I TER AT I ON 

20 CIR = 

0.93585 

DPM = 0.00026560 AT 

48 

1 1 

NSSP 

* 

ITERATION 

30 CIR a 

0.93406 

DPM a 0.00017205 AT 

48 

13 

NSSP 

X 

ITERATION 

40 CIR = 

0.93238 

DPM a 0.00019312 AT 

48 

12 

NSSP 

X 

ITERATION 

50 CIR = 

0.93085 1 

DPM a 0.00014520 AT 

47 

1 2 

NSSP 

= 


SEPARATION AT 0«1 8677 SEPARATED CP IS -0.10315 

ITERATION 60 CIR = 0.92598 DPM = 0.00040728 AT 48 12 NSSP = 


SEPARATION AT 0.15577 SEPARATEO CP IS -0.22767 

ITERATION 70 CIR a 0.91776 OPM = 0.00054908 AT 48 12 NSSP * 

SEPARATION AT 0.15577 SEPARATEO CP IS -0.19578 

ITERATION 80 CIR * 0.91135 DPM * 0.00043064 AT 47 12 NSSP * 

SEPARATION AT 0.15577 SEPARATED CP IS -0.18704 

ITERATION 90 CIR a 0.90699 DPM a 0.00042599 AT 48 12 NSSP * 
SEPARATION AT 0.15577 SEPARATED CP IS -0.18510 

ITERATION 100 CIR = 0.90266 OPM = 0.00028467 AT A7 12 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.18468 

ITERATION 110 CIR a 0.89876 OPM = 0.00027609 AT 46 12 NSSP = 

SEPARATION AT 0.1 5577 SEPARATEO CP IS -0.18447 

I TER4T I ON 120 CIR = 0.89S25 OPM a 0.00032365 AT 48 12 NSSP = 

SEPARATION AT 0.15577 SEPARATEO CP IS -0.18419 

ITERATION 130 CIR = 0.89206 OPM = 0.00022912 AT 48 12 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.18382 

ITERATION 140 CIR a 0.88916 DPM = 0.00040609 AT 48 11 NSSP = 

SEPARATION AT 0.15577 SEPARATED CP IS -0.18341 

ITER4TION 150 CIR a 0.88651 DPM a 0.00017834 AT A7 12 NSSP a 
SEPARATION AT 0.15577 SEPARATED CP IS -0.18300 

ITERATION 160 CIR * 0.88406 DPM = 0.00023943 AT 48 10 NSSP * 

SEPARATION AT 0.15577 SEPARATED CP IS -0.18259 

ITERATION 170 CIR = 0.88180 DPM a 0.00017607 AT 46 11 NSSP a 

SEPARATION AT 0.15577 SEPARATEO CP IS -0.18220 

ITERATION 180 CIR a 0.87970 DPM = 0.00025922 AT 48 11 NSSP * 

SEPARATION AT 0.15577 SEPARATED CP IS -0.18182 

ITERATION 190 CIR * 0.8777S DPM = 0.00022250 AT 48 12 NSSP = 

SEPARATION AT 0.15577 SEPARATED CP IS -0.18145 

ITERATION 200 CIR a 0.87592 OPM a 0.00014019 AT 46 12 NSSP a 
SEPARATION AT 0.15577 SEPARATEO CP IS -0.18110 

ITERATION 210 CIR » 0.67421 DPM a 0.00021975 AT 48 13 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.18078 

ITERATION 220 CIR a 0.87260 DPM a 0.00018585 AT 48 8 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.18047 

ITERATION 230 CIR a 0.87103 OPM a 0.00017825 AT 48 13 NSSP a 
SEPARATION AT 0.15577 SEPARATEO CP IS -0.18015 

ITERATION 240 CIR a 0.86964 DPM a 0.00012803 AT 46 9 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17987 

ITERATION 250 CIR a 0.86828 DPM a 0.00025523 AT 2 11 NSSP a 
SEPARATION AT 0.15577 SEPARATED CP IS -0.17960 

ITERATION 260 CIR = 0.86699 OPM a 0.00015891 AT 48 10 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17935 

ITERATION 270 CIR = 0.66576 DPM = 0.00013804 AT 48 12 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17909 

ITERATION 200 CIR = 0.86459 OPM = 0.00023842 AT 48 12 NSSP * 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17884 

ITERATION 290 CIR = 0.86348 DPM = 0.00016278 AT 48 9 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17861 


0 DELTAY QR DEL STAR a 0.1291 
0 DELTAY OR DELSTAR a 0*129! 
0 DELTAY OR DELSTAR a 0.1291 
0 DELTAY OR DELSTAR a 0.1291 
0 DELTAY QR DELSTAR a 0.1291 

0 DELTAY OR DELSTAR a 0.1428 

0 DELTAY OR DELSTAR a 0.1638 

0 DELTAY OR OELSTAR a 0.2002 

0 DELTAY OR OELSTAR a 0.2011 

0 DELTAY OR OELSTAR a 0.2041 

0 DELTAY OR DELSTAR a 0.2052 

0 DELTAY OR DELSTAR a 0.2058 

0 DELTAY OR OELSTAR a 0.2060 

0 DELTAY OR DELSTAR a 0.2061 

0 DELTAY OR OELSTAR a 0.2061 

O DELTAY QR DELSTAR a 0.2061 

0 DELTAY OR DELSTAR a 0.2060 

0 DELTAY OR DELSTAR a 0.2059 

0 DELTAY OR OELSTAR a 0.2053 

0 DELTAY OR DELSTAR a 0.2057 

O DELTAY OR DELSTAR a 0*2056 

0 DELTAY OR DELSTAR * 0.2055 

0 DELTAY OR DELSTAR = 0.2053 

0 DELTAY OR OELSTAR a 0.2052 
0 DELTAY OR DELSTAR a 0.2051 
0 DELTAY OR DELSTAR a 0.2050 
0 DELTAY OR DELSTAR a 0.2049 
O DELTAY OR OELSTAR a 0.2046 
O DELTAY OR OELSTAR a 0.20A7 



ITERATION 300 CIR = 0.862*1 DPM = 0.00013393 AT *7 12 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17837 

I TERAT I ON 310 CIR * 0.86139 OPM = 0.00010097 AT 48 13 NSSP a 

SEPARATION AT 0.’5577 SEPARATED CP IS -0.1781* 

ITERATION 320 CIR * 0.660*1 DPM = 0.00009310 AT *6 6 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17792 

ITERATION 330 CIR a 0.859*8 DPM a 0.00010669 AT *5 12 NSSP a 
SEPARATION AT 0.15577 SEPARATED CP IS -0.17771 

ITERATION 3*0 CIR a 0.85B5Q DPM = 0.000100*3 AT *7 12 NSSP * 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17752 

ITERATION 350 CIR a 0.85772 DPM = 0. 00010085 AT *6 10 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17732 

ITERATION 360 CIR a 0.85689 DPM = 0.00008178 AT *6 12 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.1771* 

ITERATION 370 CIR * 0.65609 DPM = 0.00006178 AT *8 8 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.1769* 

ITERATION 380 CIR a 0.85533 DPM * 0.00021237 AT *8 11 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17676 

ITERATION 390 CIR a 0.65*58 DPM a 0.00013*77 AT *7 8 NSSP a 

SEPARATION AT 0.15577 SEPARATED CP IS -0.17656 

ITERATION *00 CIR a 0.85387 DPM a 0.00009191 AT *6 12 NSSP a 
SEPARATION AT 0.15577 SEPARATED CP IS -0.17638 

BOUNDARY LAYER ANALYSIS FOR REYNOLDS NUMBER OF 0.63 OF 07 
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Figure 3. Separation Point and Pressure Behavior 
During Relaxation Process. 
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